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EDITORIAL 



La educacion, la cultura, la ciencia y la tecnologia son 
los elementos fundamentales en los que se basarse 
el desarrollo de un pais. Los que aplicaron estos 
elementos hace decadas el dia de hoy son las grandes 
potencias economicas del mundo. El camino ya es conocido, 
invertir en educacion y en investigacion es sentar sobre 
bases firmes nuestro desarrollo. Paises con la misma historia 
que el nuestro, como Brasil, Chile, Argentina, por citar 
algunos, estan logrando un desarrollo importante de su 
economia, por que han formulado estrategias de desarrollo 
basadas en el conocimiento, en el fomento a la innovacion 
cientifica y tecnologica para fortalecer la competitividad en 
el mercado global. Un comun denominador es que van de la 
mano el gobierno (politicas publicas), los empresarios y las 
universidades: la Triple Elice. 

Sin embargo, en Mexico no se esta actuando en este marco, 
ya que el presupuesto que se destina a las universidades 
publicas es cada vez mas escaso, y ni que decir del limitado 
apoyo que se brinda a la investigacion cientifica y tecnologica. 
El pasado mes de noviembre se realizo en la ciudad de Puebla 
la Tercera Reunion de Vinculacion de la Asociacion Nacional 
de Universidades e Instituciones de Educacion Superior del 
pais (ANUIES), en esta ocasion, participaron funcionarios 
publicos del gobierno federal, empresarios de corporativos 
importantes y universidades y tecnologicos de todo el pais. 
Lo relevante de la reunion, es que se expuso la necesidad de 
impulsar conjuntamente proyectos estrategicos, trabajar en 
redes regionales de colaboracion, formar recursos humanos 
calificados, fortalecer la innovacion y la competitividad de 
las empresas con un sentido social. En las mesas de trabajo 
se puso enfasis en el gran potencial con que cuentan las 
Instituciones de Educacion Superior en infraestructura fisica 
y recursos humanos. 

La ruta para lograr el desarrollo sustentable de nuestro 
pais ya se conoce, lo que hace falta es replicar y adaptar 
las experiencias exitosas de otros paises, integrar 
esfuerzos, destinar recursos, trabajar todos los sectores 
organizadamente. 



En este numero se presentan las cuatro secciones: 
Investigacion, Desde la academia, Politicas de ciencia y 
tecnologia y Ciencia, tecnologia y sociedad (CTS). 

En Investigacion, sobresalen los articulos que tiene que 
ver con nuestro universo, con la formacion de estrellas en 
particular, con la colaboracion de investigadores de renombre 
tanto de la Universidad de Sonora como de la Universidad 
Nacional Autonoma de Mexico y de otros paises. Se aborda 
tambien el tema de limpieza de suelos contaminados, se 
expone el diseno de un prototipo de seguidor solar, el 
problema de las caidas de adultos mayores como factor de 
riesgo, entre otros. 

La seccion Desde la academia destaca tres temas: el sismo 
de gran intensidad que se presento a principios de ano en 
el noroeste del estado de Sonora y Mexicali, la investigacion 
sobre el consumo del alcohol en los estados fronterizos y, 
finalmente, se presenta un avance sobre la planta del siglo 
XXI, la Moringa. 

El Politicas de ciencia y tecnologia se seleccionaron dos 
articulos, el primero propone unsistemadegestionambiental 
para avanzar hacia la sustentabilidad de las universidades y 
el segundo se presenta los resultados de un estudio sobre la 
gestion del conocimiento manufacturero de las pequenas y 
mediana empresas de la ciudad de Hermosillo, Sonora. 

En Ciencia, tecnologia y sociedad, se nos invita a conocer 
el posible origen astronomico de las unidades de longitud 
del sistema ingles. Un tema de interes general y de gran 
preocupacion, es el problema de la obesidad: un gran peso 
para Sonora. Un tema que nunca pierde vigencia, el agua, en 
este caso, el agua subterranea, visito desde el punto de vista 
de la fisica. 

Como es costumbre se anexan las secciones de las noticias 
mas relevantes del semestre, la gustada seccion de breviarios 
de ciencia y un reto a su conocimiento de la ciencia: El 
cientigrama. 



Una muestra del conocimiento que se genera en las 
universidades es la edicion de la revista Epistemus la cual, 
nace como una preocupacion de la Universidad de Sonora 
por divulgar el quehacer academico. Esfuerzos de este tipo 
se realizan en varias Instituciones de Educacion Superior del 
pais. 



Si usted gusta profundizar en alguno de los temas, se 
incluye el correo electronico de los autores. Escribanos, 
externe sus comentarios o sugerencias al comite editorial. 
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PARA PARTICIPAR CON ARTICULOS O COLABORACIONES EN epistemus 



La educacion, la ciencia y tecnologia son consideradas pilares 
fundamentales sobre los que se sustenta el desarrollo de un pais, 
por eso es importante fortalecer el enlace entre los que generan 
el conocimiento y los beneficiados de ello: la sociedad. Uno de los 
aspectos que distingue a la Universidad de Sonora es la generacion 
de conocimiento ysu impactoquetieneen la sociedad. Esasfcomo 
las Divisiones de Ingenieria, Ciencias Exactas y Naturales y Ciencias 
Biologicas y de la Salud crearon el proyecto editorial Epistemus el 
cual, constituye un medio de informacion y comunicacion para dar 
a conocer, en forma comprensible, las investigaciones realizadas o 
en proceso, las actividades academicas, las reflexiones en torno a 
la ciencia y tecnologia, la cultura cientifica y la sociedad, proyectos 
de vinculacion y extension de la cultura y los servicios. 

OBJETIVO 

Dar a conocer el conocimiento que se crea, genera y ensena en la 
Universidad de Sonora con finalidad de fortalecer la educacion, la 
cultura cientifica ytecnologica y de la salud de los diversos sectores 
de la sociedad. 

Areas generales de conocimiento 

• Ingenierfas: materiales, metalurgia, civil, minas, industrial, 
ambiental, hidraulica, sistemas de informacion, 
mecatronica, alimentos, energia, agua, entre otras. 

• Ciencias exactas y naturales: geologia, fisica, matematicas, 
electronica y ciencias de la computacion. 

• Ciencias biologicas y de la salud: investigacion en 
alimentos,desarrollo regional, acuacultura, salud, biologia, 
agricultura, entre otras. 

DIRIGIDAA 

A los sectores relacionados con la educacion y la investigacion, a 
empresarios, a dependencias gubernamentales, estudiantes de 
nivel medio superior y superior y sociedad en general. 

CONTENIDO DE LA REVISTA 

Articulos de proyectos de investigacion, resenas, ensayos, 
informacion sobre ciencia y tecnologia, eventos relevantes, 
convocatorias, noticias relevantes, educacion, cultura, etcetera. 

ENFOQUE DE LOS ARTICULOS 

• Los articulos reflejaran lo mas trascendente de la produc- 
cion academica, los proyectos de ciencia, tecnologia y so- 
ciedad que se realizan en la universidad y en particular en 
las tres divisiones. 

• Los articulos deberan de hacer mencion de la trascenden- 
cia de lo expuesto, su impacto en la solucion de problema- 
ticas especificas de la sociedad, del sector industrial, de la 
educacion, de la cultura, entre otras. 

• Se incluiran articulos que integren y reflexionen en torno a 
la ciencia, la tecnologia y la sociedad, que aporten elemen- 
tos precisos que permitan profundizar en el analisis y pro- 
poner esquemas de colaboracion entre los que producen 
el conocimiento y los beneficiarios o usuarios potenciales. 
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• El lenguaje escrito debe ser de divulgacion, comprensi- 
ble para un publico no especializado de nivel bachillerato 
aproximadamente, empresarios y profesionistas de otras 
especialidades. 

ARBITRAJE 

Comprende dos aspectos: 

1 . Arbitraje academico: En todos los casos, los articulos seran 
arbitrados por pares academicos. 

2. Arbitraje de divulgacion: Se procurara que el lenguaje y 
el significado de lo expuesto no pierda la idea original al 
tratarlos como articulos de divulgacion, para ello se contara 
con el apoyo en el arbitraje de academicos expertos en 
divulgacion cientifica. 

• Se les senalaran adecuaciones, modificaciones y sugerencias 
para que sean tomadas en cuenta por los autores. En todos 
los casos se definiran los tiempos para las correcciones. 

CARACTERISTICAS DE LOS ARTICULOS 

• La extension de los articulos sera maximo de 15 cuartillas, 
con figurasyfotos. 

• Incluir fotos y graficos de buena calidad en jpg o tiff de 
maxima resolucion. 

• Utilizar tipo de letra arial de 1 2 puntos con doble espacio. 

• Incluir los datos curriculares de los autores, especificando el 

area de adscripcion, perfil profesional y correo electronico. 

• Las referencias bibliograficas deberan aparecer citadas en 
el texto con un numero entre parentesis y no ser mayores 
de 3 por parrafo. 

• Al inicio del texto debera hacerse una descripcion breve del 

contenido del articulo que no sea mayor de 10 renglones 
y que de una idea clara del contenido y logre interesar al 
lector. 

• Respecto a los autores y coautores, el primero que aparezca 

sera el lider; los cuales en el caso del ensayo podran ser un 
maximo de 3 y para el caso de los resultados de investiga- 
ciones podran ser hasta 5. 

• Cuando la colaboracion sea en la modalidad de ensayo, 
hacer uso de subtitulos que faciliten la lectura del texto. 

• Para la presentacion del articulo debera de ser en dos 
versiones: electronica e impresa en papel, las cua- 
les deberan enviarse al director de la revista, al correo: 
pacheco@correom.uson.mx. 

• De publicarse el articulo, el autor debera firmar carta de 
cesion de derechos de autor, la cual se enviara al 
director de la revista una vezaprobada su publicacion por el 
comite editorial. 

FECHAS DE RECEPCION DE ARTICULOS 

• Fecha limite de recepcion de articulos: 24 de marzo de 201 1 . 

• Diseno e impresion: Junio 201 1 

Mayores informes: Con el responsable de la Division 
correspondiente. 
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MANZANARES MARTINEZ, CASTRO GARAY, ARCHULETA-GARCIA, MECTEZUMA ESQUER 



En este trabajo hacemos un resumen de la investigacion que se realiza a nivel mundial sobre las 
heteroestructuras de cristal fotonico con la idea de poner en perspectiva los resultados que se 
han obtenido sobre este tema en el Departamento de Investigacion en Fisica de la Universidad 
de Sonora. Las heteroestructuras de cristal fotonico y en general, los cristales fotonicos son una 
nueva clase de materiales artificiales con los cuales puede ser posible, en principio, desarrollar una 
revolucion tecnologica similar a la que ocurrio con la industria de la electronica a mediados del 
siglo pasado con el advenimiento del transistor. 
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CRISTALES FOTONICOS 



La explotacion tecnologica de los cristales electronicos ha 
dado lugar a una de las mas importantes revoluciones de la 
historia de la ingenieria, permitiendo un mundo moderno lleno 
de avarices tecnologicos gracias al desarrollo de la electronica. La 
mecanica cuantica que explica las propiedades de los electrones 
en diferentes materiales tambien ha permitido estudiar diversos 
aspectos de la interaccion entre la radiacion y la materia, como 
por ejemplo la emision espontanea. En el interes de controlar 
la emision espontanea que existe en ciertos materiales surgio el 
concepto de una nueva clase de materiales artificiales capaces 
de interaccionar con las ondas electromagneticas en un nivel 
estructural: los cristales fotonicos [1 ]. 

Historicamente, fue Bykov en 1972 [2] el primero en proponer 
que la emision espontanea que un cierto atomo produzca podria 
ser reducida o inhibida colocandolo en una red periodica en 
la escala de la longitud de onda de la radiacion para evitar el 
decaimiento del atomo en las bandas prohibidas. Sin embargo, 
este primer estudio no produjo mucho interes en el tema de la 
inhibicion de la emision espontanea en la comunidad cientifica. 

El verdadero despertar en el interes acerca de la inhibicion 
de la emision espontanea tuvo que esperar hasta que fueron 
publicados dos trabajos, que desde enfoques distintos 
propusieron los llamados cristales fotonicos. Estostrabajos fueron 
realizados en 1987 por Eli Yablonovitch y Sajeev John. Desde 
entonces los cristales fotonicos han sido el objeto de un intenso 
trabajo de investigacion, tanto por sus novedosas propiedades 
fisicas, como por sus potenciales aplicaciones tecnologicas [3]. 

Los cristales fotonicos usualmente estan compuestos por 
la repeticion periodica de una celda unitaria en una, dos o tres 
dimensiones en donde existen dos materiales de diferente 
indice de refraccion. En la Figura 1 presentamos una ilustracion 
de la geometria de los cristales fotonicos. En estos materiales la 
periodicidad juega un papel muy importante, ya que debido a 
fenomenos de interferencia, se espera que existan bandas de 
energia permitidas y gaps de energia prohibidos. En las bandas 
las ondas se "conducen" a traves de la red cristalina y en los 
gaps el material se considera un "aislante", ya que las ondas no 
pueden propagarse. La periodicidad de los cristales fotonicos 
generalmente esta en el rango de la longitud de onda de la luz 
visible (380-750 nm). 




Figura 1. Cristales fotonicos con periodicidad en una, dos o 
tres dimensiones. Presentamos en diferentes tonos de grises los 
segmentos con diferentes indices de refraccion. 



Los cristales fotonicos tridimensionales son analogos a las 
redes cristalinas formadas por varios elementos o compuestos 
quimicos (Si, Ge, CI, NaCI, etc) en donde la periodicidad de los 
atomos da lugar a la interferencia de los electrones con la red, 
produciendose las bandas de energia electronicas que son 
comunes en la fisica del Estado Solido. Algunos materiales 
presentan cierta forma de su estructura de bandas electronicas 
que los definen como "semiconductores electronicos" o mas 
comunmente, semiconductores. Hoy en dia la industria de la 
electronica, telecomunicaciones y la computacion esta basada 
en la explotacion tecnologica de estos semiconductores. 

Los cristales fotonicos son una nueva clase de materiales 
artificiales que pueden ser considerados "semiconductores 
opticos". Se espera que sea posible a traves de estos nuevos 
semiconductores opticos, por primera vez en la historia, lograr 
confinar, guiar y manipular el campo electromagnetico y en 
particular, la luz en forma nunca antes vista. En este momento 
estamos empezando a comprender las posibilidades que tienen 
los cristales fotonicos para la tecnologia del manana. 

En este punto es conveniente ilustrar lo importante que es hoy 
en dia el control de la luzy para ello consideremos como ejemplo 
el uso de la fibra optica para la transmision de internet. No existe 
una forma mas rapida de transmision de la informacion que la 
luz en una fibra optica. Sin embargo, aunque la informacion de 
internet se transmite por luz por medio de una fibra optica, la 
manipulacion de esta informacion se hace en las computadoras 
por medio de la electronica. El reto tecnologico de los cristales 
fotonicos es disenar dispositivos que manipulen los fotones 
en lugar de los electrones y crear, en principio, una tecnologia 
alternativa a la electronica en la que se manipulen fotones [3]. 

Para su explotacion tecnologica, los semiconductores 
electronicos han dejado una leccion: el uso de heteroestructuras. 
Una heteroestructura esta formada por dos o mas redes cristalinas 
diferentes. En la Figura 2 mostramos una heteroestructura de 
cristal fotonico. Esta heteroestructura esta compuesta por la union 
de dos redes diferentes. En el lado izquierdo, tenemos una red 
compuesta por los materiales A y B. En el lado derecho tenemos 
una red compuesta por los materiales C y D. 

Un gran ejemplo de uso tecnologico de heteroestructura 
es el transistor de estado solido, el cual dio origen a la era de la 
electronica. El transistor basico esta compuesto por tres diferentes 
redes cristalinas que se unen para lograr el control del flujo de 
electrones. Esto se logra combinando convenientemente las 
propiedades de sus diferentes bandas de energia y gracias a una 
minuscula, pero muy significativa, inyeccion externa de electrones 
para sintonizar la banda de uno de los materiales. 

En el Departamento de Investigacion en Fisica, los autores se 
han dedicado a estudiar las heteroestructuras de cristal fotonico. 
En particular, nuestra meta es la de lograr "aislantes" opticos a 
partir de heteroestructuras de cristal fotonico las cuales tengan 
una banda prohibida gigante, que no sea usual encontrar en 
cristales fotonicos canonicos. 
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A continuacion vamos a hacer una breve revision historica 
del estudio teorico de los cristales fotonicos buscando poner 
en perspectiva las contribuciones de una serie de trabajos de 
investigacion que hemos desarrollado sobre este tema. 



Cristal AB 



Cristal CD 



A B CD 

Figura 2. Heteroestructura de Cristal fotonico. En el lado 
izquierdo se muestra un cristal compuesto por los materiales 
A y B. En el lado derecho esta otro cristal, compuesto por los 
materiales Cy D. 



LA DIFICULTAD DE OBTENER CRISTALES 
FOTONICOS EN EL RANGO DEL VISIBLE 

Despues de haber sido propuesta la idea de los Cristales 
Fotonicos comenzaron los primeros estudios teoricos para 
encontrar la estructura tridimensional especifica en la que se 
pudieran encontrar efectivamente los gaps de energia prohibidos 
para ondaselectromagneticas. Los primeros resultadosteoricosse 
hicieron en forma erronea [4], ya que se utilizo una aproximacion 
escalar, muy parecida a la ecuacion de onda de Schrodinger, para 
resolver la ecuacion de onda electromagnetica y se predecian 
gaps de energia en donde en realidad no existen en el experimento. El 
problema de este primer calculo teorico radica en que cuando se busca la 
estructura de bandas en cristales fotonicos es necesario tomar en cuenta 
la naturaleza vectorial de las ecuaciones de Maxwell para plantear la 
ecuacion de onda en forma apropiada. 

La primera estructura queen forma exitosa demostro experimentalmente 
la existencia de gaps de energia prohibidas fue la Yablonovita en 1991 
[5]. La Yablonovita es una estructura cubica centrada en las caras que se 
obtiene taladrando agujeros en un bloque de material dielectrico. Esta 
estructura fue posible fabricarla en el rango de las microondas, pero no 
en el rango del visible. 

Para lograr un cristal fotonico en el rango del visible, nan sido 
propuestos diferentes metodos, como los basados en micro-fabricacion, 
autoensamblado y holografia laser. Muchas estructuras nan sido 
fabricadasen lasfrecuenciasdel microondayen el infrarrojocercano,pero 
aun sigue siendo muy dificil fabricar estructuras con bandas prohibidas 
en el visible [6]. 

iPor que es tan dificil conseguir estructuras con gap completo 
en el visible? La respuesta es extensa ya que cada estructura 
tiene sus problemas particulares, sin embargo se puede dar 
una respuesta en terminos generales; la mayoria de las tecnicas 
de fabricacion solo permiten la utilizacion de materiales de 
bajo indice de refraccion. Se ha demostrado por varios analisis 



teoricos que la estructura de bandas con un gap completo solo 
podra existir cuando se pueda construir un cristal fotonico con 
una modulacion en su indice de refraccion mas grande que 2 [7]. 
Hasta el momento no existe ninguna estructura tridimensional 
que permita fabricar estructuras con ese indice [6]. 

Los opalos artificiales alguna vez se constituyeron como una 
estructura con la cual se podria lograr estructuras en las que 
se pudiera conseguir una banda completa en el visible. Sin 
embargo, estos opalos se construyen a partir de una tecnica de 
auto ensamblado quimico, en donde unas pequenas esferas de 
material dielectrico se organizan en una estructura compacta 
cubica centrada en las caras, dejando intersticios de aire. Por 
mucho tiempo, se penso en la idea de lograr opalos inversos, 
llenando estos intersticios de aire con un semiconductor de alto 
indice de refraccion. El problema con estas estructuras es que 
existe desorden. En el 2000, Li and Zhang [8] reportaron que un 
5% de desorden estructural destruye la banda completa. Hasta 
las mejores estructuras de opalo tienen mas de esa cantidad de 
desorden. De esta forma, con los opalos nunca se ha logrado 
fabricar una estructura con gap completo en el visible 

Es necesario remarcar que el problema no es que nunca se haya 
logrado una estructura con gap en el visible. Muchas estructuras 
dan gaps en el visible. Los opalos naturales dan gaps en el visible 
y se puede observar el gap a simple vista. El problema es que 
este gap es unidireccional, en una sola direccion cristalografica. 
Usualmente la direccion TL para una estructura cubica centrada 
en las caras. El gap completo debe de ser verificado para todas 
las direcciones del espacio. Eso es un poco dificil, ya que se 
debe de cribar la estructura y lograr hacer experimentos de 
reflexion/transmision para diferentes pianos cristalinos o bien, 
experimentos con variacion de angulo de incidencia. Hasta 
ahora, ninguna estructura tridimensional ha demostrado lograr 
un gap completo en el rango del visible. 

LA INVESTIGACION ACTUAL EN CRISTALES FOTONICOS 

Aunque el estudio de los cristales fotonicos comenzo con 
estructuras tridimensionales en donde se buscaba el gap 
completo, pronto se introdujo la clasificacion de cristales 
en una (1D), dos (2D) o tres dimensiones (3D). Los cristales 
fotonicos 3D probablemente han perdido hasta el momento 
la batalla en cuanto a lograr la construccion de una estructura 
con gap completo en el visible [6]. Sin embargo, han surgido 
nuevos fenomenos de gran interes cientifico, tales como la 
supercolimacion, superprisma, refraccion negativa, entre otros 
fenomenos que estan siendo actualmente investigados en 
estructuras 3D con la idea de desarrollar nuevos dispositivos. 
Cabe decir en este punto, que el campo de los cristales fotonicos 
ha dado lugar, de alguna forma, al nacimiento o impulso de 
nuevos campos, tales como los metamateriales y la plasmonica 
[9]. 

Los cristales fotonicos en 2D son estructuras que tienen 
periodicidad en un piano y permanecen invariantes en la 
direccion perpendicular al piano. Estos cristales rapidamente 
mostraron gaps completos en el piano de la periodicidad [10] y 
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recibieron un gran impulso con Thomas Krauze en 1996 [1 1], ya 
que propuso que la mayoria de los dispositivos que se esperaba 
desarrollar a traves de los cristales fotonicos se podian hacer 
en una geometria 2D. La idea puede explicarse parcialmente 
si se considera que la tecnologia de las tabletas electronicas 
que contienen varios circuitos integrados es una "geometria 
bidimensional". Se propuso entonces que los dispositivos 
basados en cristal fotonicos solo necesitarian cristales en 2D y 
dejar el confinamiento en la direccion perpendicular al piano 
a estructuras 1D. La idea tuvo un gran impacto. La mayoria de 
la investigacion en cristales fotonicos que se hace actualmente 
es en estructuras 2D, especialmente aquellas investigaciones 
relacionadas al desarrollo de circuitos y dispositivos fotonicos. 

Los cristales fotonicos en 1 D por su parte, ya existian antes de 
que se les llamara de esa forma. Antes se les llamaba multicapas. 
Se utilizaban principalmente para hacer espejos, en especial, 
estructurasdecuartodeonda.Enlibrosclasicos,talescomoelBorn 
[1 2] se presentan metodos de calculo de reflexion y transmision 
de estas estructuras, en particular, el Metodo de Matriz de 
Transferencia. A pesar de su larga existencia, las estructuras en 
1 D siguen siendo importantes en el campo de cristales fotonicos 
ya que son estructuras sencillas en donde pueden explorarse 
fenomenos complicados. Por ejemplo, la emision espontanea 
[13], los modos de superficie [14] y los estados Tamm.Tambien 
puede considerarse que cuando se buscaba entender las 
estructuras metamateriales se comenzo con estructuras en 1D 
[15]. En nuestra investigacion sobre heteroestructuras hemos 
procedido de la misma manera y hemos desarrollado un amplio 
estudio de sistemas en 1 D, ya que apenas se esta explorando las 
posibilidades de estas configuraciones. 

Existe un punto importante que deseamos mencionar. El uso 
tecnologico de las heteroestructuas electronicas -el transistor- 
fue posible gracias a que se logro una manipulacion activa 
mediante un parametro externo de la estructura de bandas 
de uno de los materiales. En terminos simples, la corriente de 
electrones que pasa a traves de la heteroestructura puede ser 
controlada a partir de una corriente de electrones externa. Para el 
caso de los cristales fotonicos, no queda muy claro como va a ser 
posible manipular mediante un parametro externo la estructura 
de bandas. Si el parametro externo que modifica la estructura 
de bandas es la temperatura [16] o una inyeccion de portadores 
[17], estos mecanismos son muy lentos para el desarrollo de 
dispositivos fotonicos activos. Incluso la sintonizacion con 
cristales liquidos puede considerarse como un mecanismo lento 
[1 8]. Lo que se necesita es la manipulacion de luz con luz y hasta 
este momento, esa manipulacion solo parece posible a traves de 
fenomenos no-lineales. 



EL ESTUDIO DE HETEROESTRUCTURAS 

Las primeras heteroestructuras fotonicas tridimensionales 
fueron basadas en opalos artificiales [18] y diversas tecnicas de 
micro-maquinado. Desde entonces ha existido un buen numero 
de trabajos experimentales acerca del tema, aunque existen 
pocostrabajosteoricos[1 9]. Alparecer, las tecnicas defabricacion 
de opalos artificiales se adapta muy bien a la produccion de 
heteroestrucuras. Sin embargo, la heteroestructura de opalo 
hereda los inconvenientes de los opalos artificiales: el desorden 
[8]. Dado que el gap aprovechable del opalo solo existe en una 
direccion, no se ve muy claro en que forma las heteroestructuras 
de opalo puedan ser superiores a las heteroestructras en 1D. 
Por otra parte, las heteroestructuras obtenidas por tecnicas de 
micro maquinado parecen ser muy dificiles de construir [20]. 
Una estructura de micro maquinado, debido a las limitaciones 
tecnicas, usualmente es fabricada con pocos periodos. Las 
estructuras de bajo indice y pocos periodos dan como resultado 
pobres propiedades opticas. 

En la practica se exige la construccion de un gran numero de 
celdas unitarias. La investigacion sobre heteroestructuras en 2D, 
con la idea de obtener gaps gigantes es escasa. Sin embargo, 
recientemente se han reportado estructuras bidimensionales 
de bajo indice de refraccion fabricadas con tecnicas de silicio 
poroso [21] en donde seria pertinente realizar estudios para 
buscar ensanchamiento de la banda prohibida similares a los 
que hemos realizado en los calculos en estructuras 1 D. 

Las heteroestructuras unidimensionales han sido mas 
extensamente estudiadas. El primer trabajo que reporta un 
analisis sobre heteroestructuras unidimensionales fue realizado 
por H. Miyazaki et al en 1996 [22] analizando las posibilidades 
de localizacion de campo electromagnetico. El primer trabajo en 
donde se busco un gap gigante en una heteroestructura fotonica 
fue realizado en 1 998 por J. Zi et al en 1 998 [23] En este trabajo 
se estudio la superposicion de los gaps de dos redes diferentes, 
pero ambas de alto contraste dielectrico. En el 2000, estos mismos 
autores extendieron el estudio de heteroestructuras al caso de 
multiples redes cristalinas [24]. En ese mismo ano aparecieron 
los primeros resultados experimentales que reportaban la 
fabricacion de heteroestructuras fotonicas, en particular, 
fabricadas con la tecnica de silicio poroso. Desde entonces, se ha 
desarrollado una gran cantidad de articulos de investigacion en 
el campo de heteroestructras fotonicas unidimensionales, entre 
los cuales se encuentran los que se describen a continuacion. 



LAS APORTACIONES DE NUESTRO GRUPO 

Nuestra investigacion ha sido apoyada por tres proyectos del 
Consejo Nacional de Ciencia y Tecnologia (CONACYT), se han 
realizado tres tesis de Maestria en estos temas y actualmente 
esta por finalizar una tesis de doctorado. Se han publicado mas 
de 1 articulos en este tema. 
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Algunos resultados destacados se enumeran a continuacion: 



a discriminar la existencia de un gap gigante. 



• En la Referenda [25] se presenta una forma detallada 
para calcular la estructura de bandas en un cristal 
fotonico. 

• En la Referenda [26] presentamos una forma de calcular 
la velocidad de grupo de un cristal fotonico por medio de 
un calculo directo de la estructura de banda. 

• En la Referencia [27] presentamos un metodo que 
desarrollamos para estudiar heteroestructuras que estan 
construidas con materiales de alto indice de refraccion 
basado en el calculo de una supercelda, en donde, debido 
a la complejidad de la estructura de bandas, buscamos la 
forma en la que la velocidad de grupo pudiera ayudarnos 



• En la Referencia [34] presentamos un metodo para crear 
un gap gigante para heteroestructuras de bajo indice de 
refraccion. 

• En la Referencia [28] se considera la optimizacion 
de un gap omnidireccional tambien para el caso de 
heteroestructuras de bajo indice de refraccion. 

Actualmente estamos trabajando para disenar una 
heteroestructura con un gap gigante para ondas elasticas y de 
sonido. Existieron ademas algunos trabajos que no los hemos 
incluido en este resumen [29] [30] [39], pero que exploran otros 
aspectos interesantes de heteroestructuras fotonicas. 
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FORMACION DE BURBUJAS Y CASCARONES ALREDEDOR DE 

ESTRELLAS DE ALTA MASA 



BRENDA PEREZ RENDON, ALFREDO SANTILLAN, GUILLERMO GARCIA SEGURA Y LILIANAHERNANDEZ CERVANTES 




as burbujas estelares son estructuras gaseosas de gran belleza y complejidad que rodean a estrellas que 

tienen una masa mucho mayor que la del Sol. Dichas estructuras estelares son producto de una combinacion 

del medio que se encuentra entre las estrellas y del material que expulsa la estrella misma en forma de viento 

estelar. En este articulo describimos de forma basica, la fisica que hace posible la formacion de estas estructuras, 

asi como la necesidad que actualmente tiene la astronomia moderna de recurrir a simulaciones numericas para 

estudiar detalladamente la formacion de estas nebulosas. 






. 



DRA. BRENDA PEREZ RENDON. 
Departamento de Investigacion en Fisica, Universidad de Sonora. 

brenda@cajeme.cifus.uson.mx 

DR. ALFREDO SANTILLAN. 
Direccion General de Computo y Tecnologfas de la Informacion y Comunicacion, UNAM. 

alfredo@astro.unam.mx 

DR. GUILLERMO GARCIA SEGURA, 

Instituto de Astronomia, UNAM, Sede Ensenada. 

ggs@astrosen. unarm. mx 

M.I. LILIANA HERNANDEZ CERVANTES, 

Instituto de Astronomia, UNAM, Ciudad Universitaria. 

liliana@astro.unam.mx 






10 EPISTEMUS 



INTRODUCCION 

El espacio entre las estrellas contiene significantes cantidades 
de materia muy diluida, compuesta principalmente de gas 
(atomos, moleculas, iones y electrones), polvo, rayos cosmicos y 
campos magneticos, el cual esta diseminado entre las estrellas de 
forma inhomogenea, a menudo mostrando marcados contrastes 
en densidad y temperatura. Este conjunto de elementos es 
conocido como medio interestelar (MIE) y puede manifestarse en 
diferentes longitudes de onda tal como lo muestra la Figura 1. 
Para los estandares terrestres, el MIE es extremadamente tenue: 
en la vecindad del Sol, en promedio la densidad en masa es 
de aproximadamente un atomo de hidrogeno por centimetro 
cubico, lo cual es menor por 20 ordenes de magnitud que el valor 
mas bajo de la atmosfera terrestre (1). Lejos de considerar estas 
densidades tan bajas como "vacio" este gas es el responsable de 
una gran cantidad de procesos fisicos en nuestra galaxia que van 
desde la formacion de nuevas estrellas hasta la acumulacion de 
material en intrincadas estructuras dinamicas, algunas de las 
cuales son facilmente visibles con la ayuda de los telescopios. 
Estas estructuras gaseosas generalmente son esculpidas y 
alimentadas por los poderosos vientos estelares de estrellas de 
alta masa, y brillan como consecuencia de su interaccion con la 
radiacion proveniente de las estrellas. 

LAS ESTRELLAS MASIVAS Y SUS VIENTOS 

El viento estelar es un flujo de particulas que es eyectado 
por la atmosfera de las estrellas de forma continua. Es un 
mecanismo muy eficiente por el cual las estrellas pierden masa 
e inyectan energia mecanica en el medio que las rodea. Nuestro 
Sol pierde masa en una tasa de aproximadamente 10-14 masas 
solares por ano, lo que implica que si esta tasa fuera constante 
durante toda su vida, el Sol tardaria 100 billones de anos (1014 
anos) en evaporarse como consecuencia de la perdida de masa 
de su viento. La velocidad del viento solar es en promedio de 
500 km/s en la vecindad de la Tierra, alcanzando una densidad 
promedio menor a 10-30 g/cm 3 (2). En este flujo de viento solar 
esta embebido nuestro planeta y el resto del sistema solar. 
Sin embargo, la perdida de masa del Sol es muy modesta si se 
compara con los poderosos vientos estelares de estrellas de 
mayor masa, que son los responsables de moldear de forma 
dramatica las regiones que las rodean. En comparacion con el Sol, 
estas estrellas tienen tasas de perdidas de masa que son cientos 
de miles de veces mayores y sus velocidades son de cuatro a cinco 
veces mas grandes. Debido a esto, las estrellas masivas ejercen 
un fuerte impacto en la estructura y evolucion del universo a 
traves de la huella que imprimen en el MIE que las rodea. Su corto 
tiempo de vida provoca que afecten "rapidamente" su entorno, 
pues a diferencia del Sol las estrellas de mucha masa viven pocos 
millones de anos y son la principal fuente de supernovas que 
inyectan gran cantidad de energia al medio circumestelar (MCE), 
junto con elementos quimicos procesados en el interior de las 
estrellas antes y durante la explosion, en uno de los eventos mas 
extremos que se dan en el universo. 




Definimos una estrella masiva como aquella que al formarse 
tiene una masa mayor a 9 masas solares (3). Su evolucion esta 
determinada por la cantidad de material y composicion quimica 
iniciales, donde ambas condiciones establecen cuanto del 
hidrogeno que dispone puede ser utilizado como combustible 
para producir reacciones termonucleares, modificando la 
composicion quimica de la estrella. Las reaccionestermonucleares 
producen energia en el interior de la estrella que es radiada 
hacia el espacio exterior, en forma de luz y calor. En las estrellas 
masivas la nucleosintesis estelar (quemado de elementos) se 
da en el nucleo a traves de un proceso ciclico donde cada vez 
se van produciendo elementos mas pesados a partir de otros 
ligeros, que serviran de combustible para la siguiente etapa de 
quemado. Las estrellas con masas mayores a 11 masas solares 
pasan por seis etapas de quemado nuclear sucesivas en su 
nucleo: Hidrogeno ►Helio ► Carbono^Neon^Oxigeno 
► Silicio de forma no-degenerada (4), cada una de las cuales se 
enciende a una temperatura y densidad caracteristicas. Durante 
la secuencia principal, una de las etapas mas estables y de mayor 
duracion, aproximadamente 90% de la vida total de una estrella, 
se transforma en su interior hidrogeno en helio. Cabe senalar que 
el Sol se encuentra actualmente en esta etapa de evolucion. 



BURBUJAS ESTELARES 

En la secuencia principal las estrellas masivas pierden una 
gran cantidad de masa a traves de sus poderosos vientos, que 
impactan directamente sobre el MCE modificandolo quimica y 
dinamicamente. El comportamiento de esta interaccion depende 
de diversos factores como la velocidad, temperatura y masa 
del viento, las inhomogeneidades locales y/o el gradiente de 
densidad del MCE, el flujo de fotones ionizantes de la estrella, 
la tasa de enfriamiento del gas, los procesos de transferencia de 
calor y evaporacion de masa, la presencia de campos magneticos, 
etc. En algunos casos, la energia depositada en el MIE por los 
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Figura 2. La formacion de una burbuja y un cascaron de MIE barrido por un viento estelar. Izquierda. La zona [a] es donde el viento de la estrella se expande 
libremente. En [b] tenemos una burbuja de viento chocado de alta temperatura y desacelerado, que se expande por el aumento de presion termica. La region 
[c] corresponde al cascaron formado por gas barrido del MIE. I y II son los frentes de choque en reversa y directo, respectivamente. La figura fue adapatada de 
(1 3). Abajo: Grafica unidimensional de la densidad contra posicion (log p vs. r) al final de secuencia principal en la estrella de muy alta masa (30 masas solares), 
obtenida con modelos numericos nuestros. Pueden identificarse cada una de las zonas esquematizadas en la figura arriba. 






vientos estelares da lugar a la formacion de burbujas estelares. 
En la Figura 3 puede verse un ejemplo de este efecto; la Nebulosa 
de la Burbuja que se trata precisamente de una estructura creada 
por la interaction de los vientos estelares de estrellas masivas con 
el medio que las rodea. 

Realizar un estudio completo que incluya todas las variables de 
los vientos hace que el problema se vuelva muy complejo y que 
no pueda resolverse de forma analitica. Para conocer con mayor 
detalle la dinamica de la interaccion del viento con el MCE es 
necesario recurrir a simulaciones numericas que frecuentemente 
utilizan poderosas herramientas de computo. Los primeros 
intentos, antes de que hubiera computadoras suficientemente 
capaces para resolver el problema en toda su complejidad, 
fueron hechos de forma analitica proponiendo un escenario 
muy simplificado que solo arroja soluciones aproximadas al 
problema real. No obstante, estas soluciones analiticas son muy 
utiles a la hora de comprender la fisica basica de la interaccion, 
asi como, para confrontar y entender los resultados numericos. 
Por esta razon, antes de presentar nuestros resultados numericos 
discutiremos brevemente las soluciones analiticas encontradas 
por otros autores. En primera aproximacion se propone un 
escenario idealizado: un viento estelar fuerte y homogeneo que 
interactua con un medio de densidad uniforme y constante que 
rodea a la estrella, el cual como ya hemos mencionado tiene muy 
baja densidad, del orden de 10-100 part/cm3. Estos primeros 
escenarios simplificados fueron explorados por primera vez 
en las soluciones analiticas de Pikel'ner (1968) que obtienen 
una primera solucion ideal. En su modelo, un poderoso viento 
estelar barre una cavidad en el medio interestelar que la rodea. 
La configuracion resultante de este proceso esta esquematizada 
en la Figura 2 (izquierda) que muestra la configuracion del medio 
circunestelar que da lugar a la formacion de una burbuja. A la 
derecha de la misma Figura 2 se muestran los resultados de nuestro 
modelado numerico. En la cavidad barrida por el viento estelar este 
se expande libremente (region [a] de la Figura 2) y el viento finaliza 
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en un choque terminal que lo desacelera y produce una burbuja 
de viento chocado y caliente (region [b]). La presion termica de 
la burbuja caliente barre el MCE que la rodea y forma un cascaron 
denso de gas interestelar (zona [c] de la Figura 2). Otro tipo de 
soluciones analiticas autosimilares tambien fueron obtenidas por 
Weaver et al. (5) quienes resolvieron la dinamica del MCE durante 
las etapas adiabaticas iniciales, las etapas intermedias radiativas y las 
avanzadas en la evolucion, donde incorporan conduccion de calor y 
la evaporacion de material de una region a otra y obtienen soluciones 
unidimensionales que no involucran inestabilidades dinamicas. En 
sus soluciones pueden identificarse tres etapas "ideales" (5) que se 
detallaran a continuation. 

ETAPAS DE LA FORMACION DE CASCARONES Y 
BURBUJAS ESTELARES 

1) Etapa de expansion adiabatica. El viento inicial que 
emite la estrella en secuencia principal interacciona 
directamente con un medio interestelar de densidad 
constante. La evolucion es tan rapida que las perdidas 
radiativas son despreciables y el sistema se expande 
practicamente sin perdida de energia. El viento libre de 
la estrella tiene un perfil de densidad que decae con la 
distancia a la estrella (p ~ r - 2) y que acaba en un choque 
terminal, representado por un salto en la densidad que 
separa el viento libre de la estrella (region [a]) de una 
burbuja de gas caliente chocado (region [b]) que Mega a 
alcanzar temperaturas del orden de millones de grados 
Kelvin. Dentro de esta burbuja de gas caliente la presion 
es casi uniforme, como lo son tambien los valores de 
velocidad y densidad. Estas son las llamadas "burbujas 
estelares", regiones de gas chocado por el viento estelar, 
a muy altas temperaturas. Finalmente, la burbuja de gas 
chocado se encuentra rodeada por un cascaron de gas 
interestelar barrido, cuya posicion externa se expande 
en el medio interestelar. La energia mecanica inyectada 
por el viento se distribuye de diferente forma en cada 



Figura 3. Nebulosa de la Burbuja, en la constelacion de Casiopea. La figura muestra una burbuja de ( 

por los intensos vientos de la estrella central (BD+60°2522). Cabe sehalar que el cascaron que rodea a la estrella brilla en luz 

proveniente de la estrella central, asi como la region fotoionizada donde esta embebida la burbuja. Copyright: © Croman. 



bible debido a la radiacion 



region, donde 5/11 partes del total corresponden a la 
region [b], casi toda en forma de energia termica, ya que, 
la energia cinetica es practicamente nula. La otra parte se 
haya contenida en el cascaron externo: 40 % en forma de 
energia cinetica y 60% en forma de energia termica. Esta 
etapa es de muy corta duracion, tipicamente del orden 
de miles anos. 

2) Etapa radiativa. En la segunda etapa, las perdidas 
radiativas de la region [c] se vuelven comparables a las 
de la energia inyectada por el viento y el proceso de 
enfriamiento causa que el cascaron se vuelva visible en 
el optico. Las condiciones para que el material de un 
cascaron sea detectable es que alcance una densidad 
mayor a 1.67 X 10-24 gr/cm3. La masa del cascaron 
externo esta determinada por la cantidad de gas barrido. 
La Nebulosa de la Burbuja (Figura 3), situada hacia la 
constelacion de Casiopea, se trata precisamente de una 
burbuja estelar formada por este proceso. 

3) Colapso de la burbuja. En la tercera etapa, la region [b] 
empieza a enfriarse radiativamente y la burbuja tambien 
colapsa dejando que el viento de la estrella interaccione 
directamente con el cascaron externo. Sin embargo, los 
calculos en (5) indican que durante la mayor parte de la 
vida de la estrella la luminosidad radiativa es menor que 
la mecanica por lo que comunmente esta ultima etapa no 
se observa alrededor de las estrellas. 

MODELANDO NUMERICAMENTE LA FORMACION DE 
CASCARONES Y BURBUJAS 

Las soluciones autosimilares de (5) no incluyen muchas 
situaciones fisicas presentes en los cascarones reales. Por 
ejemplo, la mayoria de las estrellas masivas de tipo espectral o 
generalmente ocurren en cumulos de estrellas que estan muy 
cercanas unas de otras y el MCE de una es fuertemente afectado 
por el de sus vecinas. bien, hay estrellas en movimiento 
supersonico con respecto al medio que las rodea, que producen 
un choque de proa que modifica la estructura del cascaron 



estelar (6). Ademas de estos, hay otros procesos fisicos que no 
pueden tratarse solo con soluciones analiticas autosimilares. 

Al igual que sucedio en el campo de la evolucion estelar, 
con la aparicion de procesadores matematicos cada vez mas 
potentes, recientemente se ha puesto especial enfasis en la 
simulacion numerica de la interaccion del viento estelar y su 
entorno incluyendo nuevos escenarios fisicos, y calculando 
no solamente soluciones unidimensionales, sino en dos y tres 
dimensiones. Diversos autores han utilizado la simulacion 
numerica para el estudio de cascarones circumestelares (7), (8) 
que incluyen inhomogeneidades en densidad, inestabilidades, 
vientos episodicos, flujos de masa, etc. Todos los estudios 
concuerdan con que los cascarones creados durante la secuencia 
principal son estables aunque muy grandes y difusos por lo que 
practicamente no hay deteccion observacional de ellos. 

Las soluciones analiticas de Weaver et al. (5) solo describen el 
comportamiento del MCE cuando es afectado unicamente por 
un viento uniforme sin aceleracion. Sin embargo, en la mayoria 
de las estrellas masivas las caracteristicas del viento estelar son 
variables y se modifican a medida que la estrella evoluciona en 
su ciclo de vida. Asi pues, hay que considerar la velocidad del 
viento variable y en consecuencia, su luminosidad mecanica 
que tambien varia. La evolucion subsecuente que sera descrita 
a continuacion es resultado de una gran cantidad de estudios 
que involucran simulaciones numericas y que tratan de describir 
estas etapas de evolucion. 



a) Cascarones vinculados a Supergigantes Rojas. 

Generalmente, despues de secuencia principal una estrella 
masiva evoluciona hacia la etapa de supergigante roja (SGR) 
donde su radio se increments dramaticamente hasta alcanzar 
dimensiones de centenas de radios solares y su temperatura 
superficial decae drasticamente a un valor de aproximadamente 
3000 K. En las supergigantes rojas la tasa de perdida de masa es 
mucho mayor que en secuencia principal y el viento es lento y 
denso. Durante esta etapa, la perdida de masa por vientos es 
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Figura 4. Izquierda: La formacion del cascaron SGR dentro de la burbuja de secuencia principal. Las etiquetas son las mismas que en la Figura 2, donde 
ahadimos la zona [c'] correspondiente al cascaron SGR que esta separado del cascaron de secuencia principal por la region [b] de viento chocado. La distancia 
entre los dos cascarones no esta a escala. Notese que el viento de SGR interacciona directamente con el cascaron que forma. Derecha: Grafica unidimensional 
de la densidad contra posicion (log p vs. r) al final de SGR de una estrella de 30 masas solares), obtenida con modelos numericos nuestros. Pueden identificarse 
los dos cascarones y el resto de las zonas esquematizadas en la figura de la izquierda. 



tan grande, que la estrella Mega a perder hasta un 60% de su 
masa inicial. Este nuevo viento lento interactua con el medio 
esquematizado en la Figura 2. Despues de secuencia principal, 
el viento libre de la supergigante roja se topa con la burbuja de 
gas chocado que es menos densa cuando todavia esta en su fase 
adiabatica (region [b] de la Figura 2) y forma una sub-estructura 
de cascaron con el viento de supergigante roja. En la region [c'] 
de la Figura 4 puede verse una representacion esquematica de la 
formacion de este nuevo cascaron (izquierda). Como la creacion 
de este cascaron resulta de la interaccion de un viento lento 
despues de uno rapido (el de secuencia principal), no hay burbuja 
chocada de viento de SGR ya que el viento denso interacciona 
y forma directamente el nuevo cascaron, por lo que durante 
la etapa de supergigante roja de una estrella masiva tenemos 
dos cascarones densos separados por una region menos densa 
y caliente de viento chocado de secuencia principal, como se 
esquematiza en la Figura 4. El radio de este segundo cascaron mas 
interno (de ahora en adelante, cascaron SGR) esta determinado 
por la duracion de la estrella en la etapa de supergigante roja y 
por el balance entre la presion dinamica del viento de SGR y la 
presion termica del gas chocado originada durante la vida de 
la estrella en la secuencia principal, es decir, el tamano de este 
cascaron dependera fuertemente de la tasa de perdida de masa 
en SGR, de la velocidad del viento, de la densidad de la burbuja 
chocada de secuencia principal y de su temperatura. Al inicio de 
la formacion del cascaron SGR la presion dinamica del viento es 
menor que la presion termica de la burbuja, ocasionando que la 
frontera interna de la burbuja se mueva hacia dentro hasta que 
las presiones se igualan y el cascaron alcanza su posicion de 
equilibrio produciendo cascarones pequenos de supergigante 
roja, que son hidrodinamicamente estables (9). La posicion final 
del cascaron SGR tambien dependera de otros parametros como 
puede ser la conduccion termica entre cascarones, etc. 

Un ejemplo de este cascaron SGR puede ser el que rodea a la 
estrella de Betelgeuse en la constelacion de Orion. Actualmente 
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se considera que Betelgeuse es una estrella muy masiva, de 
aproximadamente 20 masas solares y se encuentra en su etapa 
de supergigante roja, con un radio que Mega a alcanzar alrededor 
de 900 radios solares. En 1975, A.P. Bernat y D.L. Lambert (10) 
observaron emision de una linea de potasio que parece provenir 
de un cascaron circumestelar que rodea a la estrella. El radio 
interno de dicho cascaron es de aproximadamente 10 radios de 
Betelgeuse. Algunasestrellasmasivas despues depasar por todos 
los ciclos de quemado nuclear en su interior producen un nucleo 
de hierro y explotan como supernovas tipo II, durante su etapa 
como supergigante roja, como presumiblemente sucedera con 
Betelgeuse. La explosion y el material estelar eyectado impactara 
directamente sobre el cascaron SGR formado previamente por la 
estrella. 

b) {Que pasa con las estrellas muy masivas? 

Cuando las estrellas son mas masivas que las descritas en 
la seccion anterior, siguen una evolucion posterior a la de 
supergigante roja, que dependera de su masa inicial, metalicidad, 
etc. Cuando las estrellas nacen con masas superiores a 20-25 
masas solares, despues de pasar por la etapa de supergigante 
roja, la estrella evoluciona nuevamente hacia un estado de 
menor radio y mayor temperatura aumentando la velocidad del 
viento (~ 500 km/s). Algunas de estas estrellas que se encuentra 
en la etapa de supergigante azul (SGA), finalmente explotan 
como supernova tal como ocurrio en el caso de 1 987A (1 1 ). Otras 
estrellas despues de la etapa de supergigante roja se transforman 
directamente en estrellas tipo Wolf-Rayet (WR), que son estrellas 
masivas (con masas mayores a 20-30 masas solares), de muy altas 
temperaturas y que se encuentran en etapas muy avanzadas 
de su evolucion. Se caracterizan principalmente por sus fuertes 
vientos estelares a traves de los cuales pierden una gran cantidad 
de masa. Dependiendo de el subtipo, muestran intensas lineas 
en emision de hidrogeno y helio ionizado, e incluso carbono, 
nitrogeno y/o oxigeno. En la etapa WR la estrella se despoja de 
su envolvente de hidrogeno a traves de un viento muy masivo 




y rapido (del orden de 2000 km/s) que impacta sobre el MCE 
que ha sido modificado por las etapas de secuencia principal y 
supergigante roja. 

Independientemente de que la estrella se convierta en SGA o 
WR, se genera nuevamente un viento rapido que forma un nuevo 
cascaron dentro del viento de lento de SGR. Durante esta tercera 
etapa se forma un tercer cascaron (que llamaremos cascaron azul 
cuando es formado por el viento de SGA y cascaron WR cuando 
es formado por el viento de WR). La formacion y dinamica de 
estos cascarones es muy compleja y ha sido tratada de forma 
numerica por nuestro equipo de trabajo (12). La velocidad del 
tercer cascaron es significativamente mas grande que la de 
viento y que la del cascaron de SGR, por lo que se expande mas 
rapidamente y colisiona con este ultimo. Las circunstancias y 
consecuencias de esta colision pueden ser encontradas en (9). 
En la Figura 5 se muestran 2 ejemplos de burbujas estelares, NGC 
6888 y la S308, que presentan cascarones de gas alrededor de 
estrellasWolf-Rayet. 



Por otro lado, en la Figura 6 presentamos un ejemplo de los 
resultados de una de nuestras simulaciones numericas vinculados a 
la formacion de varios cascarones en 2D, que son parte de nuestros 
esfuerzos por entender la cinematica y dinamica alrededor de estos 
cascarones. Pueden observarse dos estructuras de burbuja/cascaron 
alrededor de una estrella Wolf-Rayet, situada en el centro de la 
figura. Estas estructuras son objetos celestes relativamente escasos 
debido a la corta duracion de su vida y al hecho de que provienen 
de estrellas masivas que no son las mas comunes en nuestra galaxia. 
Sin embargo, su observacion nos provee de una oportunidad 
unica de estudiar estos raros objetos, lo cual incluye confrontar las 
caracteristicas observadas con el modelado numerico que se realiza 
para explicar su naturaleza. Este campo de estudio en la astronomia 
es actualmente un campo muy dinamico, que cambia conforme 
vamos adquiriendo y desarrollando mayores herramientas de 
computo que han resultado indispensables para el estudio de este 
y muchos otros fenomenos astronomicos. 
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Figura 6. Simulacion numerica de un cascaron formado alrededor de una 
estrella Wolf-Rayet. La imagen muestra en niveles de gris el logaritmo de la 
densidad (g/cm 3 ) del material que rodea a la estrella antes de que estalle como 
supernova. Puede observarse la posicion del cascaron de WR a una distancia de 
aproximadamente 1 pc (32.6 ahos luz), el cual resulta fragmentado debido a la 
accion de inestabilidades dinamicas en el gas circunestelar. 
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INTRODUCCION 

En las diversas areas de la ingenieria se presentan ciertos tipos 
de problemas que podriamos llamar arquetipos ya que sirven 
para ilustrar conceptos, tecnicas y procedimientos en una forma 
particularmente clara e integral. Uno de estos arquetipos en 
ingenieria electrica/mecanica lo es el problema del control de un 
pendulo invertido. Este tipo de pendulo, que tiene una variedad 
de presentaciones, se le puede encontrar en diversos sistemas 
de la vida cotidiana. Uno de los mas conocidos es el relativo al 
juego de mantener la vertical de una escoba invertida apoyada 
sobre la palma de la mano (figura 1). Es bien sabido que esto se 
logra mediante los desplazamientos rapidos y oportunos de la 
posicion de la mano en la direccion por donde se aprecia que la 
escoba se inclina. Esta habilidad solo se consigue con un poco 
de practica y solo en el caso de que la escoba no sea muy corta. 




Figura 1 Equilibrando una escoba 

A nivel laboratorio, el sistema de pendulo invertido mas 
sencillo que se puede construir consiste de una varilla anclada 
sobre una superficie en un punto de apoyo sobre el cual la varilla 
puede girar en un piano. La superficie o movil puede desplazarse 
en la direccion coplanar al piano de giro de la varilla. Para lograr la 
vertical, al movil se le debe de aplicar una fuerza con la magnitud 
y direccion adecuada que contrarreste la tendencia de la varilla a 
caer.Vease la figura 2. 

En las siguientes lineas trataremos de explicar primero en 
que consiste el problema del pendulo invertido y sus aspectos 
centrales. Enseguida se analiza la dinamica del sistema, es decir, 
aplicamos las leyes de la mecanica para determinar con precision 
como es el comportamiento del sistema bajo ciertas condiciones. 
Despues se explica con cierto detalle como se implementa un 




Figura 2 Pendulo invertido sobre movil y la fuerza F 

controlador en base a la informacion de la dinamica del sistema. 
Posteriormente se describe como se construyo un modelo real. 
Finalmente se abunda un poco mas sobre diversos aspectos del 
controlador del pendulo invertido. 



EL PROBLEMA DEL CONTROL 

Puesto que implementar un control automatico para equilibrar 
una escoba equivale a desarrollar un control para un pendulo 
invertido, nuestro enfoque es analizar y construir un modelo de 
pendulo invertido en la forma de una varilla giratoria tal como la 
que aparece en la figura 2. 

Para desarrollar un mecanismo automatico que permita 
controlar esta varilla (pendulo) se requieren tres cosas esenciales: 

• Una forma continua de medir el angulo del pendulo 
respecto a la vertical. 

• Un movil impulsado por un mecanismo con una fuerza 
variable. 

• Un procedimiento que relacione, de alguna forma, la fuerza 
del mecanismo impulsor con la posicion y velocidad del 
pendulo. Este es un procedimiento de retroalimentacion ya 
que la variable de salida, en este caso el angulo, modifica a la 
variable de entrada, la fuerza impulsora. 

Intuitivamente pareceria claro que con desplazar el movil en la 
direccion opuesta al angulo del pendulo este podria mantener el 
equilibrio. Esto lo hacemos, mas o menos, al controlar la escoba 
invertida. Sin embargo, en un sistema de control automatico 
tenemos que determinar con toda precision la cantidad de 
impulso que el sistema tiene que aplicar al movil en todo 
instante. Se debe disponer de una formula o relacion que le 
permita al sistema determinar la cantidad del impulso aplicado. 
Cuando equilibramos la escoba, el cerebro "calcula", de algun 
modo, la cantidad de impulso que debemos aplicar en la mano 
para lograr el equilibrio. Esta cantidad no es constante pues varia 
dependiendo como se"siente"la escoba. El exito no se logra sino 
despues de cierto entrenamiento que le permitan al cerebro, 
nervios y musculos probar diferentes tecnicas de equilibrio. 
Ciertamente no podriamos describir con precision cual es este 
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algoritmo que el cerebro emplea. Es algo que a la conciencia le 
queda muy difuso. 

Existen varias tecnicas para determinar la regla o algoritmo 
que un sistema automatico puede seguir para controlar con 
exito un sistema. El metodo clasico, que enseguida mostramos, 
requiere conocer en detalle la dinamica del sistema. Esto es, 
partiendo de las caracteristicas del sistema (masa, dimensiones, 
etc.) y de las leyes de la fisica que lo rigen, tratar de predecir el 
comportamiento del sistema ante determinados estimulos. 
Teniendo esta informacion podemos intentar alguna estrategia 
cuantitativa- o sea, que se pueda expresar matematicamente- y 
evaluarla teoricamente para anticipar su resultado. Por lo tanto, 
lo primero que se tiene que averiguar es como varian la posicion 
y velocidad del pendulo y del carro cuando este ultimo se suelta 
a partir de un angulo inicial y se somete a una cierta fuerza F. 

DINAMICA DEL SISTEMA 

En la figura 3 se muestra el esquema de fuerzas que actuan 
sobre el sistema y que permitiran estudiar su dinamica. 



mgL<p — mL —^- = (J + mL 2 ) — ? 




Figura 3 Parametros y esquema de fuerzas del sistema 

Las variables y parametros del sistema son: 

X: Posicion horizontal. 

Y: Posicion vertical 

0: Angulo de la varilla respecto a la vertical 

2L: Longitud total de la varilla 

M: Masa del movil (carro) 

m: Masa de la varilla 

F: Fuerza aplicada al movil 

La descripcion de la posicion y velocidad de la varilla y el 
movil se obtiene al aplicar las leyes de Newton al sistema. Las 
ecuaciones que se generan considerando pequenas desviaciones 
angulares [1] son: 
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En donde J=mL 2 /3 es el llamado momento de inercia de la 
varilla y g es la gravedad. Estas ecuaciones diferenciales lineales 
describen como se mueven el carro y el pendulo si conocemos 
sus condiciones iniciales, es decir, la posicion y velocidad del 
pendulo y el carro en un momento inicial dado. Se obtienen al 
sumar las componentes horizontales y verticales de las fuerzas 
y las torcas que actuan en el sistema. Son ecuaciones lineales de 
segundo orden cuya solucion se puede obtener por metodos 
clasicos o cuya simulacion se puede hacer en computadora por 
medio de programas como MATLAB. Por supuesto, la naturaleza 
de la fuerza aplicada F es crucial en la forma de respuesta que 
tendra el sistema. Y es aqui justamente en donde tenemos que 
aplicar la idea del control automatico. Como se menciono al 
principio, la idea basica para poder controlar el pendulo (o la 
escoba) consiste en aplicar una fuerza en la direccion contraria a 
la direccion en que va cayendo el pendulo. En nuestro problema 
el carro solo se puede desplazar en la direccion del eje X y, por 
la forma en que se articula el pendulo, este solo puede rotar en 
el piano X-Y. Esto simplifica la accion de control porque para 
mantener la posicion vertical del pendulo solo tenemos que 
desplazar el carro en la direccion del eje X: en sentido positivo o 
en sentido negativo. El unico problema es determinar que tanto 
se debe desplazar el carro para lograr el objetivo del control. Y 
como la posicion del carro la determina la fuerza F aplicada, el 
problema se resume en la pregunta: {De que magnitud debe ser 
F para lograr el control del pendulo? 



EL CONTROLADOR 

Existen varias respuestas a esta pregunta pero la mas simple, y 
quizas obvia, es la de hacer que la magnitud de F sea proporcional 
al angulo del pendulo pero con signo contrario. Es decir 



F = -k* 



Esto es lo se conoce con el nombre de control proporcional. 
La siguiente pregunta obligada es: <;Esto es suficiente para 
lograr la verticalidad del pendulo? La unica forma de responder 
adecuadamente a esta pregunta es analizar matematicamente 
el sistema incorporando la ecuacion (3) a las ecuaciones (1) y (2) 
y resolver para O(t). Antes de seguir procede una observacion: 
si la posicion inicial de la varilla es perfectamente vertical (O=0) 
entonces la fuerza aplicada sera cero y la varilla permanecera 
indefinidamente equilibrada. Esto es, por supuesto, en la teoria. 
En la practica sabemos que cualquier pequeno disturbio la sacara 
de su posicion. Por lo tanto, para ser mas realista, se supondra 
que inicialmente la posicion de la varilla es diferente de 0=0. De 
ahi en adelante el sistema buscara la posicion de equilibrio. 

No entraremos en los detalles de como resolver las ecuaciones 
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Figura 4 Respuestas de la varilla al control proporcional 

1, 2, y 3. Solo mencionaremos que en la teoria del control se 
acostumbra emplear laTransformada de Laplace [2] para analizar 
y resolver estas ecuaciones. En particular, con esta tecnica se 
desarrolla el concepto de Funcion de Transferencia [3] que 
permite visualizar los rasgos generales del com porta miento del 
sistema. Para este caso particular se puede demostrar [1] que se 
tienen tres tipos de respuestas diferentes dependiendo del valor 
particular de la constante k. Para el rango de k < g(M+m) la fuerza 
es demasiado debil para sostener el equilibrio de la varilla y esta 
cae. Para el rango de k > g(M+m) la fuerza es demasiado grande 
y la varilla oscila alrededor de la vertical. Solo en el caso de que 
k sea exactamente igual a g(M+m) la varilla se mantendra en la 
posicion inicial del angulo. De hecho, si la posicion original no 
es cero, esta se mantendria indefinidamente. En la grafica de la 
figura 4 se observan las curvas de respuesta correspondientes a 
estos tres casos. 

El tercer caso en el que k = g(M+m) es impracticable. Cualquier 
desviacion, por pequena que sea, nos llevara a alguno de los 
otros casos. Por lo tanto, podemos concluir que para una fuerza 
como la de la ecuacion (3) el sistema es inestable y no se puede 
lograr la verticalidad de la varilla. Puede quedar la impresion 
de que aumentando el valor de k la varilla pudiera alcanzar la 
vertical. Pero intuitivamente podemos ver porque no es asi: para 
grandes valores de k la recuperacion de la varilla sera tan rapida 
que se nos ira en la direccion contraria; para cuando cambie de 
signo la fuerza, la varilla reaccionara en direccion contraria pero 
con demasiada rapidez. El resultado neto es que la varilla oscilara 
rapidamente alrededor de la vertical un cierto angulo que no 
disminuira con el tiempo. 

MEJORANDO EL CONTROLADOR 

En el anterior intento por controlar el pendulo nemos pasado 
por alto un detalle que no es obvio pero que tiene logica: nos 
estamos preocupando unicamente por aplicar una fuerza de 
recuperacion que solo depende del tamano del angulo 0, que 
es el error entre el angulo deseado (0 grados) y el angulo actual- 
no estamos tomando en cuenta la velocidad angular con la que 
va cayendo la varilla, o sea dflVdt. Por ejemplo, si la varilla esta 
desviada de la vertical un angulo pequeno pero la velocidad 
con la que esta cayendo es grande, la fuerza F sera pequena y 



no anticipara que pronto aumentara la desviacion. Un instante 
posterior el error se habra incrementado y ya no recuperara 
lo ganado. Podemos mejorar la informacion sobre la fuerza 
modificando la relacion (3) para que incluya un termino de 
velocidad. La nueva propuesta quedaria asi: 

F=<krp+k 2 f) 

En esta nueva relacion la magnitud de la fuerza F depende 
tanto de la desviacion (D como de la velocidad angular 60/ 
dt. Los factores de peso son k1 y k2 que se conocen como 
constante proporcional y constante derivativa respectivamente 
y a la relacion (4) se le conoce como controlador proporcional 
+ derivativo o simplemente controlador PD. Queda la cuestion 
de que valores les podemos dar a estas constantes. En la teoria 
del control lo que se analiza es la forma que en que se modifica 
la funcion de transferencia del sistema con el agregado del 
controlador PD a fin de escoger los mejores valores de k1 y k2 que 
satisfagan las especificaciones requeridas. Equivalentemente, se 
puede incorporar la ecuacion (4) con las ecuaciones (1) y (2) y 
buscar lo valores de las constantes que produzcan un sistema 
estable, si es que los hay, por supuesto. Afortunadamente las 
ecuaciones se pueden resolver por metodos clasicos ya que son 
lineales y de segundo orden. Al hacerlo, se puede demostrar [3] 
que la solucion para produce un sistema estable si se cumplen 
las condiciones k2 > y k1 > Mg. 
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Figura 5 Respuesta del sistema con control PD 

En la figura 5 se aprecia como es la respuesta del sistema para 
valores de k1 = 250 y diferentes valores de k2. 

Se aprecia que el sistema sigue produciendo oscilaciones pero, 
a diferencia del sistema proporcional, ahora las oscilaciones se 
amortiguan conforma pasa el tiempo. Entre mayor es el valor de 
k2 mas rapido se amortiguan. Esto es, la varilla tiende a la vertical 
mas rapidamente. 
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UN MODELO REAL 

En la practica, la construccion de un pendulo invertido y su 
correspondiente sistema de control presenta una gran variedad 
de dificultades. La principal tiene que ver con la forma en que 
se genera la fuerza F con las caracteristicas de la ecuacion (4). 
En el Laboratorio de Electronica del Departamento de Fisica 
hemos empleado un carrito electrico de juguete impulsado por 
un pequeno motor con su transmision integrada para construir 
el modelo del pendulo invertido. Desmantelamos la carroceria 
dejando solo el chasis como plataforma para soportar la varilla/ 
pendulo. En su base la varilla se articula al eje de un sensor de 
posicion tipo incremental que se encuentro sujeto al chasis del 
carro y que sirve para medir el angulo $. Vease la figura 6. 

Sobre el mismo chasis se localiza el circuito que se diseno 
para efectuar la accion de control PD. Este circuito, basado en un 
microcontrolador, lee constantemente el sensor de posicion de 
la varilla y asi determina (|) y d(|)/dt en todo momento. Con estos 
valores y empleando un algoritmo para el calculo de la ecuacion 
(4) determina el voltaje que debe aplicar al motor del carro. Este 
voltaje esta directamente relacionado con la fuerza F por lo que 
la accion de control se completa. 

El modelo construido presenta el inconveniente del 
deslizamiento o patinaje de las ruedas sobre el piso; esto 
degrada la accion de control. Pudimos resolver el problema, 
hasta cierto punto, utilizando una cubierta de plastico rugoso 
sobre el piso formando una pista por donde el carro efectua 
sus desplazamientos. Otras aproximaciones mas eficientes en 
este aspecto emplean rieles con bandas de traccion para mover 
el carro [4], [5]. Obviamente, esto introduce complicaciones 
mecanicas y limita el espacio de operacion. 

MAS SOBRE EL CONTROL 

Una pregunta natural es si existe otra expresion diferente de 
la ecuacion (4) que permita obtener resultados satisfactorios, 
es decir, que produzca un sistema estable. La respuesta es que 
hay muchos. Nada mas que el control PD es, para este sistema, 
el mas sencillo de todos. Una de las modificaciones que se le 
pueden hacer a la expresion (4) es agregarle un termino que 
incluya la integral del angulo (|) con una cierta constante de 
proporcionalidad. Esto da origen a un controlador muy empleado 
en sistemas de control llamado control PID. De hecho, nuestro 
controlador es de este tipo pero con la constante de la integral 
igual a cero. Normalmente este termino sirve para disminuir el 
error de la variable de salida cuando se alcanza el estado estable. 

En principio, cualquier cambio es factible hacer al controlador 
y estudiar sus efectos en la estabilidad del sistema. El problema 
es que cuando las ecuaciones diferenciales que caracterizan al 
sistema son de mayor orden, el analisis del sus comportamiento 
es mas complicado y su implementacion fisica tambien lo es. Por 
ejemplo, en nuestro caso del pendulo invertido, si incorporamos 
el efecto del motor del carro en el sistema, a las ecuaciones (1) 
y (2) habra que agregarles otra ecuacion diferencial de segundo 
orden que incluye las variables de torque y corriente electrica; 
con esto el sistema se convierte en uno de orden cuarto que 
complica bastante su analisis. No se debe de olvidar que estas 
ecuaciones son validas para valores pequenos de (|), del orden de 
200 o menos. En caso contrario, las ecuaciones se convierten en 
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Figura 6 El carro con el pendulo invertido 

no lineales y las tecnicas de control mencionadas no funcionan 
adecuadamente. 

Existen diferentes aproximaciones al problema del control de 
sistemas cuando estos son muy complicados o se desconocen sus 
ecuaciones caracteristicas. Una de las mas exitosas es la conocida 
con el nombre de control difuso [6] que se acostumbra emplear 
cuando se tiene informacion empirica sobre el comportamiento 
del sistema sin necesidad de conocer el detalle de su estructura 
interna. 



CONCLUSIONES 

El problema de mantener equilibrada una escoba con la palma 
de la mano ilustra muy bien el problema tecnico del control de 
un pendulo invertido. Se han mostrado los detalles de como 
implementar un control automatico para este sistema que 
permiten comprender mejor los conceptos de retroalimentacion, 
control proporcional y control derivativo. Se ha mostrado como 
es indispensable conocer la dinamica del sistema para escoger los 
parametros correctos del controlador. La implementacion fisica de 
un pendulo invertido utilizando un movil en una dimension (como 
el desarrollado en nuestro laboratorio) es un clasico experimento 
que se puede emplear en un curso formal de teoria de control. Los 
detalles tecnicos de su construccion y los calculos involucrados 
aportan elementos didacticos de gran valor. 



BIBLIOGRAFIA 

1 . Elgerd Olle I, Control Systems Theory, McGraw-Hill ISE, Tokyo, 1 967. 

2. Spiegel Murray R, Ecuaciones Diferenciales Aplicadas, 3a Edicion, 
Prentice-Hall, Mexico, 1 983. 

3. Ogata Katsuhiko, Modern Control Engineering, Prentice Hall, New 
Jersey 1997. 

4. Stang John, The inverted pendulum, Cornell University Engineering 
Division, Design Project Report, May 2005 

5. Medrano-Cersa, G.A, Robust computer control of an inverted 
pendulum, IEEE Control Systems Magazine, Vol. 1 9, No. 3, Junio 1 999. 

6. Wang Li-Xin, A Course In Fuzzy Systems and Control, Prentice-Hall, 1 996 



INVESTIGACION 




DISENO Y ELABORACION DE PROTOTIPO SEGUIDOR 

SOLAR DE DOS GRADOS DE LIBERTAD PARA CELDAS 

SOLARES FOTOVOLTAICAS 



GUSTAVO JESUS OZUNA HUERTA, MARIA ELENA ANAYA PEREZ, NUN PITALUA DIAZ, 

CARLOS ANAYA EREDIAS 



El presente articulo propone y evalua mediante simulacion matematica, un sistema difuso de seguimiento 
solar con el objetivo de lograr el correcto posicionamiento de una celda solar fotovoltaica y asi obtener la 
mayor cantidad de luz solar durante el dia y por lo tanto producir la mayor cantidad de energia electrica. 
Los resultados de la primera fase de la investigacion muestran el comportamiento de los motores de corriente 
directa con respecto a los sensores propuestos para desarrollar esta funcion. Esta propuesta es una alternativa 
en la mediana generacion de energia electrica para satisfacer la demanda de electricidad en casa habitacion o 
comunidades ru rales de forma sustentable. 



ING. GUSTAVO JESUS OZUNA HUERTA 
gozuna@industrial.uson.mx 

M.C MARIA ELENA ANAYA PEREZ 
meanaya@industrial.uson.mx 

DR. NUN PITALUA DIAZ 
meanaya@industrial.uson.mx 

M.C CARLOS ANAYA EREDIAS 

Departamento de Ingenieria Industrial 

Universidad de Sonora 



EPISTEMUS 21 



INTRODUCCION 



PROTOTIPO 



El presente articulo es un producto del proyecto de la 
investigacion "Diseno y Elaboracion de Prototipo Seguidor Solar 
de Dos Grados de Libertad Para Celdas Solares Fotovoltaicas"en 
el cual se propone un mecanismo para eficientar la captacion de 
energia solar mediante el posicionamiento de una celda solar 
con respecto al sol, el cual en su primera fase se modelo el control 
de este dispositivo matematicamente con algoritmos de control 
difuso, que permiten garantizar el correcto posicionamiento 
del panel fotovoltaico durante el dia, y asegurar la maxima 
captacion solar para producir mayor cantidad de energia 
electrica. Los algoritmos de control posicionan el panel de este 
a oeste siguiendo la trayectoria aparente del sol y de norte a sur 
posicionandolo en el cambio de estaciones en el ano, como se 
muestra en la Figura 1. 





w 

s 




^/"v 4 ^^--/ 



Figura 1. Posiciones del sol y del panel solar 



El objetivo de la investigacion es mostrar una alternativa mas 
en el seguimiento solar desarrollada totalmente en el campus 
y mostrar la viabilidad de esta por simulacion matematica. La 
estructura del siguiente articulo muestra la problematica a 
resolver , seguida del diseno del prototipo, el modelado del 
motor de corriente continua que utilizara el prototipo, el sistema 
de control y los resultados de la simulacion, asi como la propuesta 
electronica por parte del equipo de investigacion que da pie a la 
segunda fase de este proyecto de investigacion. 



PROBLEMATICA 

Hoy en dia podemos observar a varias disciplinas que buscan 
desarrollar sistemas sustentables de transformacion de energia 
y que abordan el problema de varios puntos de vista para lograr 
una produccion mas limpia. Las celdas solares fotovoltaicas que 
podemos observar en la ciudad estan instaladas fijas y captan 
en su mayoria los rayos del sol de medio dia. Ademas en el 
mercado hay diferentes alternativas de seguimiento solar para 
celdas fotovoltaicas como el sistema de posicionamiento global, 
utilizando controladores logicos programables y mediante 
computadora utilizando LabView, el prototipo propuesto busca 
ser una alternativa de bajo costo. 



En el sistema mostrado en la Figura 2 se necesita posicionar 
ambos motores, para asegurar el correcto seguimiento del panel 
solar en cualquier posicion del sol con respecto al Este-Oeste o 
Norte-Sur. Los 2 motores que se utilizan son de corriente directa, 
con la siguiente aplicacion: 




Figura 2. Seguidor solar con sensores 



Ml: Motor 1 — > Posicion horizontal (Norte-Sur, Sur-Norte) 

M2: Motor 2 — ^Posicion vertical (Este-Oeste, Oeste-Este) 

Estos son modelados matematicamente con las siguientes 
ecuaciones: [4], [5] 









Donde Vapp es el voltaje aplicado, R es la resistencia electrica, 
L es la inductancia, b es la friccion, J es el momento de inercia 
para el motor y KO es la contante del campo electromagnetico 
representados en la Figura 3. 




Figura 3. Motor DC 
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Ademas para el sistema se estan considerando 5 sensores 
codificados de la siguiente manera: 

51 -> Este 

52 -> Oeste 

53 -* Sur 

54 ~> Norte 

55 -> DiaoNoche 

Estos sensores foto resistivos solo permiten pasar un de- 
terminado voltaje que va desde los volt hasta los 5volts 
dependiendo de la cantidad de luz solar que les llegue, logrando 
con ello valores continuos. El quinto sensor S5, solo es sensible a 
la radiacion solar e indica si es de noche o de dia. 



CONTROL 

Para el control se propone utilizar un controlador difuso debido 
a que este emula aceptablemente el razonamiento humano 
y podra tomar decisiones a partir de datos imprecisos [2]. Se 
pretende utilizar un Microcontrolador en el cual se programen 
las reglas difusas y asi poder posicionar el panel fotovoltaico en 
los cuatro puntos cardinales como lo muestra la Figura 1 

Las ecuaciones del motor DC asi como las entradas de voltaje 
por parte de los sensores fueron simuladas por medio de 
MATLAB y Simulink lo cual es representado en la Figura 4.Para 
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Figura 4. Sistema de control de posicion para el adecuado seguimiento 

del sol en Simulink 



el sistema de control se utilizan los dispositivos de entrada S1 y 
S2 para el sistema de control de posicion del motor 1, S3 y S4 
para el sistema de control de posicion del motor 2 y ademas el 
sensor S5 el cual detecta la radiacion solar y controla los otros 
cuatro sensores. El Microcontrolador (simulado) realiza las 
comparaciones de los sensores S1- S2 y S3-S4 para calcular el 
error y despues enviar la serial al motor 1 o motor 2, lo cual es 
retroalimentado constantemente, mediante la comparacion del 
voltaje de entrada de los sensores para lograr el posicionamiento 
del panel fotovoltaico, adecuado si el error es cero cuando los 
sensores S1 y S2 esten recibiendo el mismo voltaje. 

Los sensores S3 y S4 trabajan de la misma manera. La posicion 
se determina con la informacion establecida en las tablas 
de inferencia difusa (tablas 1 y 2), donde cada regla hace la 



comparacion entre el sensor S5 y el error horizontal (eH) para 
el caso del motor 1 y el error vertical (eV) para el motor 2. 
Las funciones de membrecia utilizadas para el sistema de 
seguimiento solar difuso estan representadas en las Tablas 1 y 2: 



eH\S5 


ON 


OFF 


NEG 


EO 


OK 


Ok 


OK 


OK 


POS OE 


OK 



Tabla I.Tabla de inferencia difusa para posicion horizontal 



ev\S5 


ON 


OFF 


NEG 


SN 


OK 


Ok 


OK 


OK 


POS 


NS 


OK 



Tabla 2. Tabla de inferencia difusa para posicion vertical 




Figura 5. Conjunto difuso de entrada para ambos errores de posicion (eH, eV) 

En la Figura 5 podemos observar la representacion del 
conjunto difuso del voltaje resultante en eH, donde a cada valor 
de voltaje de [-5 5], le corresponde un rango [0 1] segun sea su 
pertenencia. Cuando el sistema queda estabilizado en "OK", el 
sistema de seguimiento recibe en su entrada la misma cantidad 
de voltaje entre los sensores y por lo tanto la celda fotovoltaica 
esta posicionada correctamente en la trayectoria Este-Oeste. 
Ademas el conjunto difuso para eV es similar al utilizado en la 
figura 5 porque el comportamiento de los errores son parecidos, 
debido a que los voltajes de entrada en ambos sistemas de 
control provienen de sensores con el mismo rango de voltaje. 




Figura 6. Conjunto difuso de salida para el 
controlador difuso del motorl 
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Figura 7. Conjunto difuso de salida del sensor 5 



Figura 9. Trayectoria linguistica para la posicion horizontal del panel solar 



Los conjuntos difusos de salida que se muestran en las figuras 
6 y 7, son similares para los motores 1 y 2, puesto que ambos 
comparten las mismas caracteristicas mecanicas y electricas 
ademas de recibir la informacion por parte de sensores con el 
mismo com porta miento 

En la figura 7 se muestra el conjunto difuso de entrada del 
sensor S5 el cual es un pirheliometro y su funcion es detectar 
la radiacion solar. Este se ubica fuera del sistema, por lo cual se 
podrian controlar varios sistemas de seguimiento solar difuso a 
la vez. Cuando es de dia S5 toma el valor ON, por lo que los otros 
sensores empezaran a buscar la maxima luz solar y cuando el 
sensor S5 esta en OFF, los otros sensores dejan de buscar, aunque 
los sensores del panel detecten otro tipo de luz, por ejemplo de 
lamparas u otras fuentes luminosas, no seran tomadas en cuenta 
porque solo es sensible a la radiacion solar, evitando con esto 
perdidas de energia durante la noche estando en reposo. 

En las tablas 3 y 4 se muestran las reglas de control difuso para 
el sistema tanto en posicion vertical y horizontal. 



FU ZZY RULE S 
1 ' ir ieH is NEG i and ;£D rs ON i Ihen (Ml is EW) 



2 If («H Is ok * and (85 fe ON) Than |M1 Is OK> 



"1 flNeH is PQ5) and (55 Is OM> then (Ml b WE) 



if icH is NEG) and i S5 is OFF) ttwn {Ml is OK) 



5 If (eH is ok.) and fSS is OFF) (hm (M 1 is OK) 

6 " IMeH ra POS) and (EA is UF-hHhen (Ml is UK) 



Tabla 3. Reglas base para el control difuso de 
posicion horizontal 





FUZ2Y R ULES 


7 


IT (ev b MEG) and (3d is ON) inert ( M2 13 sm> 


B 


If (cV is ok) and (S5 is ON? own (M2 is OK) 


i 


IT<eV is F'OS) ana |S& is ON) tnEn (MZ is Nt>) 


ib 


1 HeV is NL<j J airnl (S6 q OF F ) thon (MZ * UK) 


11 


it (ev & ok) and (35 \& or r ) in#n \ m? >s ok* 


it 


ITfeV ft FOS) and (55 ft OFr)tn«l iM2 ft OK) 



Tabla 4. Reglas base para el controlador difuso de 
posicion vertical 
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SIMULACION 

Para la simulacion se representaron los valores de voltaje por 
medio de 5 vectores como entrada para los sensores S1, S2, S3, 
S4 y S5 en el sistema representado en la Figura 4 por medio del 
software SIMULINK perteneciente a la plataforma MATHLAB. La 
Figura 9 muestra el comportamiento del motor 2 que controla 
la posicion de la celda solar, en la zona WE-EW (Oeste-Este Este- 
Oeste), donde la posicion del panel comienza en la zona WE y 
luego va en direccion de EW para direccionarse hacia el punto 
maximo de EW para despues regresar a la zona WE terminando 
en el punto maximo de WE. 
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Figura 10. Velocidad angular horizontal en un periodode 12 hrs 

Asi mismo, la figura 10 muestra el comportamiento del 
motor 1, y la posicion del panel durante el dia y en la zona NS- 
SN (Norte-Sur Sur-Norte). La grafica comienza en la zona NS el 
comportamiento de este vector siempre se mantiene en el rango 
de [5 - 5] lo cual indica que la trayectoria se va intercambiando 
de las zonas NS-SN. La cual nos muestra la trayectoria linguistica 
por medio de zonas de activacion en un periodo de 12 hrs. Las 
Figuras 10 y 11 muestran la velocidad angular horizontal del 
motor 1 y la velocidad angular vertical del motor 2 en un periodo 
de12hrs[1]. 



PROPUESTA ELECTRONICA 

Despues de haber realizado las pruebas de simulacion del 
comportamiento de los motores con respecto a los sensores y 
haber obtenido que el sistema simulado sea estable, el equipo 
de investigacion elaboro la propuesta electronica, la cual se 
muestra en la figura 12. Este circuito utiliza el Microcontrolador 




Figura 1 1 . Velocidad angular vertical en un periodo de 1 2 hrs 
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Figura 12. Propuesta electronica para el seguidor solar 



PIC 16F628A el cual tendra programadas las reglas difusas 
ademas de ser un Microcontrolador popular de bajo costo. 

CONCLUSIONES 

El buscar nuevas alternativas para eficientar los sistemas 
tradicionales de transformacion de energia solar en la universidad 



de Sonora, contribuira al logro de una universidad sustentable 
comprometida con el medio ambiente. El sistema propuesto es 
una alternativa economica en la mediana generacion de energia 
electrica para satisfacer la demanda de electricidad en aulas, casa 
habitacion o comunidades rurales. 
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CAIDAS EN EL ADULTO MAYOR 

FACTORES DE RIESGO QUE REPERCUTEN EN SU ESTADO DE SALUD 



BLANCA ESCOBAR CASTELLANOS, MARIA GPE. MARTINEZ MARTINEZ, MARIBEL VILLEGAS CASTRO 



Las caidas en los adultos mayores (AM) son un problema de salud publica porque se encuentran 
vinculadas a la discapacidad, a la disminucion de la calidad de vida y al incremento del numero de 
anos de vida potencialmente perdidos. Ademas, las lesiones demandan una mayor inversion de 
recursos economicos, sociales y humanos, la planeacion de programas especializados y la redefinicion 
de prioridades en salud. 

El presente trabajo permitio identificar los factores de riesgo que contribuyen a que un adulto mayor 
presente una o mas caidas en el curso de su vida, se identifico que los mas susceptibles de presentarlas 
eran los mayores de 75 anos. Entre los factores intrinsecos se identifico la perdida del equilibrio como 
principal causa. Una minoria mostro riesgo a la marcha identificandose mecanismos de adaptacion 
para la misma, como el uso de aditamentos y aumento de la base de sustentacion para deambular. En la 
valoracion del equilibrio-marcha como predictor de caidas se encontro que un porcentaje importante 
presentaba cinco veces mayor riesgo de caer. 

Entre los factores de riesgo extrinsecos, las causas de caidas mas frecuentes fueron atribuidas a la 
presencia de obstaculos y superficies resbaladizas dentro del entorno habitual de la poblacion, sin 
embargo las caidas se presentaron dentro y fuera del hogar. Ademas se identifico que la mayoria de la 
poblacion no cuenta con la cultura de prevencion de caidas. 






MCE. BLANCA ESCOBAR CASTELLANOS. 

Departamento de Enfermeria. Universidad de Sonora. 

Correo:escobarcblanca@enfermeria.uson.mx 

M. ED. MARIA GPE MARTINEZ MARTINEZ. Facultad de Enfermeria. 
Universidad Autonoma de Nuevo Leon. Correo: mmartinez20@hotmail.co 

MCE. MARIBEL VILLEGAS CASTRO. 

Departamento de Enfermeria. Universidad de Sonora. 

Correo: maribelvillegas_72@hotmail.com 



26 EPISTEMUS 



INTRODUCCION 



Entre los accidentes, las caidas son la fuente mas importante 
de morbilidad y mortalidad para los Adultos Mayores (AM), 
convirtiendose en un grave problema de salud publica, por las 
lesiones que producen. De los ingresados en el hospital a causa 
de las caidas, unicamente el 50 % estara vivo al cabo de un ano 
(1). En paises desarrollados el 30 % de las personas de 65 anos o 
mas que viven en su casa, pueden caer una o mas veces por ano; 
este porcentaje se incrementa hasta 50% en los adultos mayores 
de 80 anos. El fenomeno de traumatismos originados por caidas 
predomina en el sexo femenino en una relacion 2.7:1.00 en 
las edades de 60 a 65 anos, pero tiende a igualarse por sexos 
conforme se incrementa la edad y tiende a ser mas frecuente en 
los varones despues de los 85 anos (2). 

Segun la Encuesta Nacional de Salud en Mexico (2000) dentro 
de las causas de los accidentes, las caidas ocuparon cerca de la 
mitad de los casos notificados. La prevalencia en las mujeres se 
incremento al menos dos puntos porcentuales por cada grupo 
decenal de edad a partir de los 60 anos, pasando de 4.6 % en el 
grupo de 60 a 69 anos, a 6.6 % en el grupo de 70 a 79 anos, hasta 
alcanzar 8.9 % en el grupo de 80 y mas anos (3). 

La caida es la consecuencia de cualquier acontecimiento que 
precipita al individuo al suelo en contra de su voluntad, siendo la 
etiologia multifactorial, por su elevada frecuencia y sus multiples 
complicaciones constituyen un importante problema de salud. 
Son la principal causa de muerte por lesion en los AM debido a la 
susceptibilidad a los accidentes, atribuido a cambios fisiologicos 
o enfermedades degenerativas agudas o cronicas (4) 

La interaccion entre factores de riesgo (alteraciones del 
equilibrio y marcha, farmacos y enfermedades cronicas) y la 
exposicion a situaciones u oportunidades de caer se deben tomar 
en cuenta simultaneamente, debido a que la caida en un adulto 
mayor no es un hecho impredecible o inevitable, debido al azar 
o al envejecimiento normal del individuo. En la mayoria de los 
casos se trata de una inadaptacion entre la persona y su entorno 
de origen multifactorial; el riesgo de sufrir una caida aumenta 
linealmente con el numero de factores de riesgo existentes (5). 



CAUSAS Y CONSECUENCIAS QUE PRECIPITAN AL 
ADULTO MAYOR A PRESENTAR UNA CAIDA 

Algunos autores identificaron que las caidas en la poblacion 
de AM, son causadas como muestra de lo que en apariencia 
pareceria ser el comienzo de un proceso de claudicacion de la 
fuerza vital de un sujeto, asociando el proceso de envejecimiento 
con varios cambios que predisponen a las caidas: los responsables 
de la estabilidad postural, las alteraciones visuales, auditivas 
y neurologicas propias de la edad (6). En relacion al lugar de 
ocurrencia aproximadamente 80 % ocurren en el hogar y el resto 
fuera de el; suceden tanto en posicion de pie como en situacion 
de reposo (7). De un tercio a la mitad de las personas mayores de 
65 anos sufren una caida cada ano (8). 

Aproximadamente 30 % de los adultos mayores que se 
encuentran en casa hogar para ancianos sufren caida cada ano; 
el numero es mayor en las instituciones de salud. A pesar de que 
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menos de una cafda de cada 10 causa una fractura, un quinto de 
los incidentes de cafdas requiere atencion medica (9). Es comun 
que los AM que presentan cafdas tengan en promedio de tres 
a cuatro enfermedades coexistentes, entre las mas comunes 
se encuentran alteraciones cardiovasculares, metabolicas, 
musculo esqueleticas, asf como discapacidad visual y motora 
predominantemente en el sexo femenino y de mayor edad (2,6). 
Las cafdas repercuten en la estancia hospitalaria y elevan el costo 
economico por las complicaciones inherentes (10). 



FACTORES DE RIESGO 

Para este estudio se consideraron los factores de riesgo tanto 
intrinseco como extrinseco que pueden propiciar una cafda y 
afectar la integridad fisica de las personas mayores de 60 anos. 

Los factores de riesgo son cualquier caracteristica o 
circunstancia detectable de una persona o grupo de personas 
que se sabe asociada con un aumento en la probabilidad de 
padecer, desarrollar o estar especialmente expuesto a un proceso 
morbido. Estos factores de riesgo (biologicos, ambientales, 
de comportamiento, socioculturales, economicos) pueden 
sumandose unos a otros, aumentar el efecto aislado de cada uno 
de ellos produciendo un fenomeno de interaccion (11,1 2). 

El factor de riesgo intrinseco se considero como una condicion 
fisiologica relacionada con el individuo que puede predisponer 
a una caida, como son: cafdas previas, numero de enfermedades 
y medicamentos que actualmente consumen las personas asf 
como el nivel cognitivo, equilibrio y cambios en la marcha. 

El factor de riesgo extrfnseco, se reconocio como la condicion 
ambiental que puede incrementar la susceptibilidad de caer y 
producir danos a la salud del adulto mayor, como: caracterfsticas 
de los pisos, de las areas de la vivienda, la iluminacion, escaleras, 
areas exteriores y precauciones personales. 



OBJETIVOS DEL ESTUDIO 

El proposito de este estudio fue identificar los factores de riesgo 
de cafdas a los cuales esta expuesta la poblacion adulto mayor 
derechohabiente de una institucion de primer nivel de atencion. 
Los objetivos fueron: 1. Identificar el perfil sociodemografico y 
antecedentes de cafdas de la poblacion sujeto de estudio y 2. 



Identificar los factores de riesgos intrfnsecos y extrfnsecos de 
cafdas de la poblacion por sexo y grupo etareo. 

CARACTERISTICAS DEL ESTUDIO Y METODOLOGIA 

El presente estudio es de tipo descriptivo, a traves del cual 
se describen los factores de riesgo de cafdas de los AM (13). 
La poblacion estuvo conformada por adultos mayores de 60 
anos adscritos a la Clfnica de Medicina Familiar del ISSSTE de la 
ciudad de Hermosillo, Sonora, Mexico. La muestra fue calculada 
en el paquete estadfstico nQuery advisor v. 4.0, con un nivel 
de significancia de 0.05 y un Ifmite de error de 0.06, quedando 
constituida por 190 AM, correspondiendo a 107 hombres y 83 
mujeres. 

Para recolectar la informacion, se acudio al servicio de 
consulta externa de la institucion de salud para la seleccion de 
los usuarios con base al control de citas diarias. Se elaboro una 
cedula de datos sociodemograficos la cual permitio recabar 
informacion sobre: edad, sexo, escolaridad, problemas de salud, 
tratamiento medico actual, asf como antecedentes de cafdas, 
causas y lugar de ocurrencia. Se valoro el estado cognitivo con 
el examen mental de Pfeiffer. El equilibrio estatico y dinamico fue 
evaluado con la Escala de Tinetti para equilibrio. Continuando 
con la valoracion de la marcha empleando la Escala de Tinetti 
para marcha. Finalmente se utilizo una lista de Verificacion de 
Seguridad Domiciliaria (14) la cual evaluo las caracterfsticas 
arquitectonicas y la presencia o ausencia de riesgos en el 
medio ambiente del adulto mayor. Al concluir la entrevista se 
proporciono informacion sobre Prevencion de Cafdas en el Hogar 
de los Adultos Mayores. 

Los resultados se procesaron en el paquete estadfstico SPSS v. 
13.0, se analizo la distribucion de variables continuas mediante 
la prueba de Kolmogorov Smirnov. Se obtuvo el coeficiente 
de Alpha de Cronbach de los instrumentos de Tinetti para las 
subescalas de equilibrio y marcha y la escala equilibrio-marcha. 
Se utilizo estadfstica descriptiva a traves de frecuencias y medidas 
de tendencia central. La estadfstica inferencial fue analizada con 
pruebas no parametricas de U de Mann Whitney y la Prueba de 
Kruskal-Wallis. 

RESULTADOS 

La edad predominante de los adultos mayores estuvo en el 
rango de 65 a 74 anos ( = 71.43 + DE = 7.08), el 56.3 % era del 
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sexo masculino y 43.7 % femenino, el 88.4 % vivia con familia 
y 85.3 % habitaban en casa de una planta. El 55.8 % habia 
presentado antecedentes de caidas previas al momento del 
estudio, con = 1.14 + (DE =1.54) de estos, 42.1 % sufrio de una 
a dos caidas y 7.9 % de cinco a mas, un 32.1 % habian sufrido 
alguna caida en los ultimos seis meses. La poblacion padecia de 
una a dos enfermedades cronico-degenerativas (70.5 %) y el 75.2 
% tomaban de uno a cuatro medicamentos ( = 2.60 + DE = 1 .77). 

De la poblacion que reporto haber sufrido caidas, el 41.6 % 
correspondio a factores extrinsecos predominando la presencia 
de obstaculos y superficies resbalosas en el medio ambiente. 
Del 14.2 % de los factores intrinsecos prevalecio la perdida del 
equilibrio. En cuanto al lugar de ocurrencia de las caidas 26.3 
% sucedieron en la casa, 16.8 % en la calle y 12.6 % en lugares 
publicos. De las caidas en casa predomino el area del bano (6.8 %) 
y patio (4.7 %). Las caidas registradas en la calle sucedieron en la 
banqueta (14.2 %). Las caidas en lugares publicos prevalecieron 
los comercios en 2.1 % y otras con 6.9 %. 

Al evaluar el estado cognitivo se identifico que 91 % de los 
adultos mayores no presentaron alteracion y solo 6.3 % tenia 
deterioro leve. En la evaluacion del equilibrio 89.5 % mostro 
riesgo de perdida del equilibrio y 44.2 % presento riesgo para 
la marcha. Con relacion a la evaluacion integrada de equilibrio - 
marcha se identifico que mas del 50 % de la poblacion presento 
algun riesgo de caer y 20.5 % presento cinco veces mas riesgo. 
(Tablal). 



Variable 


Indicador 


f 


% 


Estado Cognitivo 


Intacta 

Deterioro Leve 
Deterioro Mode ra do 
Deterioro Severo 


173 
12 

4 
1 


91.1 
6.3 
2.1 
0.5 


Equilibrio 


Normal 
Con riesgo 


20 
170 


10.5 
89.5 


Marcha 


Normal 
Con riesgo 


106 

84 


558 
44.2 


Equilibrio-Marcha 


Normal 

Riesgo 

5 veces mas riesgo 


53 
98 
39 


27.9 
51.6 
20.5 


Fuente: SPMSQ, Tinetti 




n= 190 



Tabla 1. Evaluacion del Estado Cognitivo, Equilibrio, Marcha y Equilibrio- 
Marcha 

En la valoracion del equilibrio, 66.3 % empleo los brazos para 
levantarse, 21.1 % requirio de mas de un intento o ayuda para 
ponerse de pie, en posicion de pie 68.9 % mantuvo una posicion 
estable con una distancia mayor de 12 centimetros entre los 
talones, al ejercer un empujon ligero 45.8 % se tambaleo y 14.7 
% tendio a caerse, 42.1% al mantenerse de pie con los ojos 
cerrados mostro inestabilidad, al realizar un giro de 360 grados 
con los ojos cerrados 44.2 % mostro pasos discontinuos y 33.2 % 
inestabilidad, al sentarse el 44.7 % uso los brazos para asegurarse 
y 7.4 % se mostro inseguro. 

En la evaluacion de la marcha, el 88.9 % no vacilo al iniciar la 
marcha, presentaron buena altura, longitud y simetria del paso 
con ambos pies, 1 3.7 % mostro discontinuidad al caminar, 23.6 % 
presento desviacion leve/moderada o marcada de la marcha, en 



23.7 % se observo alguna alteracion de la posicion del tronco y 
29.5 % separo los talones cuando camino. 

Con respecto a los factores de riesgo extrinseco en el interior 
del hogar el 24.7% contaba con pisos de material resbaloso, 37.4 
% presentaban desnivel o irregularidades. En el bano no utilizaba 
tapetes antideslizantes en el area de regadera (84.2%), ni en el 
cuarto de bano (70.5 %), un 86.3 % no tenian instalada agarradera 
fija en el area de regadera. En relacion con la iluminacion 45.8 % 
no contaban con luces nocturnas que conducen al bano y 51 .6 % 
con lampara o apagador que pudieran alcanzar desde la cama. 

En cuanto a los factores de riesgo extrinsecos en el area 
exterior, el 56.8 % tenian tapete en la entrada de la casa, en 
30.5 % no estaban bien pavimentadas las calles y banquetas de 
su comunidad y 22.6 % contaban con obstaculos alrededor de 
su casa. En relacion al calzado 28.4 % acostumbraba a caminar 
descalzo y con calcetines (25.3 %). Entre las precauciones 
personales el 65.3 % no tenian banco para alcanzar objetos 
en altos y 41.6 % tenian animales de compania con los cuales 
pueden tropezar. Con respecto a la poblacion que habitaba en 
casa de dos plantas 12.1 % no contaban con antiderrapante en 
cada uno de los escalones. 

Al comparar el sexo de los participantes con los factores 
de riesgo intrinsecos se identifico diferencia significativa 
con mediana mas baja para el sexo masculino en relacion a 
enfermedades (Mdn= 1.00, U= 3103.50, p= .000), tratamientos 
(Mdn= 2.00, U= 31 67.50, p= .001 ) y caidas previas (Mdn= 0.00, U= 
271 7.00, p= .000). Ademas se identifico diferencia significativa en 
el sexo femenino en relacion al equilibrio con mediana = 1 1.00 
(U= 361 7.00, p= .028). (Tabla 2). 



Factores de riesgo 
intrinsecos 




Sexo 




U 


Valor 
dep 


3 Femenino 


6 Masculino 


.V 


Mdn 


X 


Mdn 


Enfermedades 


2.00 


2.00 


1.38 


1.00 


3103.50 


.000 


Tratamientos 


3-07 


3.00 


2.23 


2.00 


3167.50 


.001 


Caidas previas 


1.71 


1.00 


0.69 


0.00 


2717.00 


.000 


Equilibrio 


10.61 


11.00 


11.50 


13.00 


3617.00 


.028 


Marcha 


10.31 


12.00 


10.32 


12.00 


4368.00 


.832 


Equflibrio-marcha 


20.96 


22.00 


21.82 


24.00 


3729.00 


.058 


Fuente: CDS, Tinetti 
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b n = 107 



Tabla 2. Factores de Riesgo Intrinsecos en relacion al Sexo de los 
participantes 

Los factores de riesgo intrinsecos (enfermedades, tratamientos 
y caidas previas) en relacion al grupo etareo, no se encontro 
diferencia significativa, entre los grupos de edad con p> .05. 
En relacion con enfermedades se encontro en el grupo de 75 
a 79 anos mediana mas alta (Mdn= 2.00, H= 1.85, p= .763), en 
tratamientos la mediana mas baja se identifico en el grupo de 
75 a 79 anos (Mdn= 2.00, H= .25, p= .992), y en caidas previas la 
mediana mas baja fue para el grupo de 65 a 69 con mediana de 
0.00 (H= 5.09, p= .277). 

En cuanto a los factores de riesgo intrinsecos (equilibrio, 
marcha y equilibrio-marcha) en relacion con el grupo etareo, 
se identifico diferencia significativa entre los diferentes grupos 
de edad con p= .000. Obteniendose medianas mas bajas en el 
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equilibrio para los grupos de 75 a 79 y los mayores de 80 anos con 
mediana = 1 0.00 (H= 29.1 1 , p= .000) respectivamente. En marcha 
el grupo de 80 y mas anos mostro mediana mas baja (Mdn= 
1 0.00, H= 1 8.04, p= .001 ) y en relacion con el equilibrio-marcha el 
grupo de 80 y mas anos mostro mediana mas baja (Mdn= 19.00 
(H=28.56,p=.000).(Tabla3). 





Factor de Riesgo 












Valor 


Intrinseco 


Edad 





X 


Mdn 


H 


dep 




60-64 


31 


13.16 


13.00 








65-69 


48 


12.63 


13.00 






Equilibrio 


70-74 


48 


10.50 


12.00 


29.11 


.000 




75-79 


33 


10.18 


10.00 








80 y mas 


30 


8.60 


10.00 








60-64 


31 


1 1 .29 


12.00 








65-69 


48 


11.27 


12.00 






Marcha 


70-74 


48 


10.25 


12.00 


18.04 


.001 




75-79 


33 


9.94 


12.00 








80 y mas 


30 


8.30 


10.00 








60-64 


31 


24.35 


25.00 








65-69 


48 


24.02 


24.50 






Equilibrio-Marcha 


70-74 


48 


20.75 


23.00 


28.56 


.000 




75-79 


33 


20.12 


22.00 








80 y mas 


30 


16.90 


19.00 






Fuente: Tinetti 
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190 



Tabla 3. Factores de Riesgo Intrinsecos: Equilibrio, Marcha y Equilibrio- 
Marcha en relacion al Grupo Etareo 

Al valorar los factores de riesgo extrinsecos en relacion con 
el sexo se identifico que no mostraron diferencia significativa, 
la mediana mas baja correspondio al sexo femenino (Mdn= 
67.18, H= .781, p= .377). Con respecto a los factores extrinsecos 
en relacion con el grupo etareo, se observa que no existio 
diferencia significativa entre los grupos de edad, las medianas 
mas bajas se observaron para los grupo de 60 a 64 y 80 y mas 
anos respectivamente (Mdn = 67.1 8, H= 5.237, p = 264). (Tabla 4) 

DISCUSION 

La mayoria de la poblacion AM estudiada vivia con familia y 
mas de la mitad habia presentado caidas en algun momento 
de su vida, de estos casi la mitad habian sufrido de una a dos 
caidas en un periodo mayor a seis meses y una tercera parte en 
los ultimos seis meses, esto coincide con autores (2) quienes 
senalan que en los paises desarrollados las personas mayores 
de 65 anos que viven en su hogar pueden caer una o mas 
veces por ano, incrementandose hasta la mitad en los AM de 80 

30 EPISTEMUS 



anos, identificando como factor de proteccion que el AM viva 
acompanado. 

Mas de la mitad de los AM de 60 anos tenian antecedentes de 
caidas previas, el predominio de las mismas ocurrieron fuera del 
hogar atribuidas a obstaculos y superficies resbalosas; algunos 
autores mencionan que la mayoria de las caidas ocurren fuera 
del hogar ocasionadas por las mismas causas, sin embargo, 
este hallazgo no coincide con lo reportado por otros autores (7) 
quienes identificaron que la mayoria de las caidas ocurren en el 
hogar, tanto en posicion de pie como en situacion de reposo. 

En relacion al lugar especifico donde ocurrieron las caidas 
dentro del hogar, los factores ambientales son la principal 
causa precipitante de caidas por accidente siendo los mas 
representatives baneras, escaleras y muebles que disminuyen 
las areas de transito en los pasillos del hogar, mismos que 
se constituyen en factores de riesgos de caidas. Los AM 
padecian hasta cuatro enfermedades cronico degenerativas, 
predominando hipertension arterial sistemica y diabetes 
mellitus y recibian en promedio 2.6 hasta un maximo de nueve 
medicamentos como tratamiento, esto es similar a lo reportado 
por algunos autores (2,6) donde es comun la presencia de 
pluripatologias, medicacion multiple, edad mayor de 80 anos, 
discapacidad para deambular y visual como factores asociados 
a caidas previas. 

En relacion al estado cognitivo, la gran mayoria de la poblacion 
tenia un nivel cognitivo intacto, no se encontro relacion con las 
caidas; estos resultados coinciden con autores que reportaron 
que las personas que cayeron tenian evaluacion funcional 
y mental normal en dos terceras partes de la poblacion, no 
encontrandose asociacion con el riesgo de caer. En la evaluacion 
equilibrio-marcha se identifico que la mitad de las personas 
mostraron riesgo y una quinta parte presento cinco veces mas 
riesgo de caer, los trastornos de marcha estan relacionados con el 
riesgo de caer atribuidos a cambios fisiologicos de la edad. 

Las personas que presentaron alguna alteracion del equilibrio 
y/o marcha, mostraron mecanismos de adaptacion para su 
situacion como son el empleo de los brazos para levantarse 
y asegurarse al momento de sentarse, algunos requirieron 
de ayuda para ponerse de pie y mantuvieron una base de 
sustentacion con distancia mayor a 1 2 centimetros al ponerse de 





Variable Edad 


11 


X 


Afdn 


H 


Valor 
dtp 


60-64 

c . 65-69 
Factores d© 7Q 74 

Riesgo Extrinsecos 

3D y mis 


31 

4B 

33 
30 


66.33 
68.39 
71.02 
69.27 
69.11 


67.18 
6875 
7031 
68.75 

67.18 


5237 


.264 


Fuente: LVSD 
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Tabla 4. Factores de Riesgo Extrinsecos en relacion al Grupo Etareo 

pie. Para la marcha se identificaron el uso de ayudas para caminar 
y separar los talones cuando caminan. Esto coincide con otros 
autores (1 , 8) que reportan que a los factores de riesgo se agregan 
el no flexionar las rodillas durante la marcha y no levantar el pie 
completamente del suelo. 

En relacion con los factores de riesgo extrinsecos se identifico 
que la mayoria de la poblacion estaba expuesta a riesgos de caer, 
esto coincide con lo reportado por algunos autores (9,14) que 
encontraron asociacion entre las condiciones de la vivienda con 
el riesgo de caer e identificaron que la mayoria de la poblacion 
mostro algun tipo de riesgo de caer dentro y fuera del hogar, 
siendo los mas detectados la presencia de pisos resbalosos, 
desniveles, obstaculos y el uso inapropiado de calzado. 

CONCLUSIONES 

La poblacion estudiada estaba expuesta a riesgo como 
antecedentes de una a dos caidas previas en un periodo de 
seis meses, siendo los mas susceptibles de presentar caidas los 
mayores de 75 anos. Entre los factores intrmsecos se identifico 
que la perdida del equilibrio fue la principal causa atribuida a las 
caidas. En razon de las enfermedades, consumo de medicamentos 




y caidas previas se encontro diferencia significativa en relacion 
con el sexo. 

En la valoracion del equilibrio se encontro que casi la totalidad 
de la poblacion presento alteracion o riesgo de caer, en la 
marcha una minoria mostro riesgo, se identificaron mecanismos 
de adaptacion para la misma, como fue el uso de aditamentos 
y aumento de la base de sustentacion para deambular. En la 
valoracion del equilibrio-marcha como predictor de caidas 
se encontro que la mayoria presentaba riesgo y un porcentaje 
importante presentaba cinco veces mayor riesgo de caer. 

En relacion con los factores de riesgo extrinsecos, las causas 
de caidas mas frecuentes fueron atribuidas a la presencia de 
obstaculos y superficies resbaladizas dentro del entorno habitual 
de la poblacion, sin embargo las caidas se presentan dentro y 
fuera del hogar. 

No se encontro diferencia significativa en los factores de riesgo 
extrinsecos en relacion con el sexo y el grupo etareo. Al realizar 
la verificacion de seguridad domiciliaria y medidas de seguridad 
personal se identifico que la mayoria de la poblacion no cuenta 
con la cultura de prevencion de caidas. 
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DINAMICA CLASICA Y FRAGMENTACION EN EL PROCESO DE 
FORMACION DE PROTO-ESTRELLAS 



GUILLERMOARREAGA GARCIA, SILVIO ORESTE TOPA 



Las estrellas nacen en nubes de gas como resultado final del proceso de colapso gravitacional. 
En este proceso actuan varias fuerzas, a saber: la prominente fuerza atractiva de la gravedad 
debido al peso de la propia nube; la fuerza repulsiva de la presion del gas y la fuerza centrifuga 
originada por la rotacion que se observa generalmente en las nubes moleculares. En este 
trabajo de divulgacion nos interesa describir algunas de las principales caracteristicas dinamicas de 
este proceso, haciendo enfasis en la ocurrencia de fragmentacion de la nube progenitora. Para esto, 
nos apoyaremos en simulaciones numericas del proceso de formacion de proto-estrellas que hemos 
realizado recientemente. 
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INTRODUCCION 



El proceso de formacion estelar ocurre en nubes de hidrogeno 
molecular que se localizan frecuentemente en los brazos de las 
galaxias espirales. En la Figura 1 se muestra una fotografia de 
la nebulosa NGC 1333, una nube progenitora tipica de estrellas 
ubicada en la constelacion de Perseus a 1000 anos luz de 
distancia, una de las regiones mas proximas de formacion estelar. 
Otro ejemplo muy conocido es la nebulosa M42, o nebulosa de 
Orion. A unos 1 270 anos luz de la Tierra, es una de las nebulosas 
mas brillantes del cielo nocturno, observable facilmente a simple 
vista al sur del cinturon de Orion en la constelacion del mismo 
nombre. Los astronomos nan concluido por observaciones que 
las estrellas recien nacidas estan agrupadas muy frecuentemente 
en pares y en menor proporcion en grupos de mas de dos 
miembros de estrellas. En el seno de estas nubes gigantes 
puede haber entonces miles de estrellas jovenes apareadas en 
pequenos grupos binarios. 

Ahora bien, podemos considerar que el proceso de formacion 
estelar viene ocurriendo desde que el Universo era muy joven 
y actualmente sigue teniendo lugar en diferentes partes del 
Cosmos, dando origen a millones de estrellas jovenes. Sin 
embargo, existen ciertas caracteristicas comunes en estos 
conjuntos de estrellas jovenes, por ejemplo: el numero de 
estrellas recien nacidas con una cierta masa parece que siempre 
es muy similar independientemente de la nube progenitora en 
la cual se realice el conteo. Esta similitud en la distribucion de 
estrellas jovenes podria ser una consecuencia de que el proceso 
de formacion de estrellas es similar en todas las nubes que se 
encuentren en condiciones fisicas similares en los brazos de 
cualquier galaxia. En otras palabras, deben existir principios 
fisicos generales que rigen el complejo proceso de evolucion que 
va desde la formacion de la nube progenitora hasta el nacimiento 
de estrellas. 




Figura 1. Nube Molecular en NGC1333 mostrandovariasgeneracionesde estrellas. 



Si tenemos en cuenta las escalas de tiempo involucradas en el 
proceso de evolucion estelar - millones a miles de millones de 
anos - la pregunta que naturalmente surge es <;c6mo podemos 
comprender el mismo si nuestra existencia es tan efimera en 
comparacion?. Haciendo uso de una analogia, podemos imaginar 
a un cientifico que solo tiene un mes para estudiar una abundante 
poblacion de leones en el Africa. Es comprensible que en ese 
periodo de tiempo le sea imposible poder observar la evolucion 
de un individuo en particular desde se gestacion hasta su muerte 
por vejez. La solucion al problema viene del hecho de que debido 
a que la poblacion es abundante, tendra ejemplos de los distintos 
estadios de estos animales (leonas prenadas, cachorros, animales 
adultos y en el final de su vida). De ello se desprende que el 
cientifico pueda extrapolar sus datos y reconstruir la vida de esta 
poblacion de animales sin la necesidad de observar la historia 
de un solo animal. Un astronomo al observar el firmamento se 
encuentra en la misma situacion, pudiendo recrear la evolucion 
de los distintos tipos de estrellas y elaborar hipotesis sobre como 
las leyes fisicas operan sobre la misma. De alii la importancia de 
las simulaciones numericas, ya que nos permiten comparar y 
validar o refutar estas hipotesis en base a los resultados de las 
mismas y los datos recopilados observacionalmente. 

El trabajo de investigacion que hemos desarrollado en anos 
recientes nos ha permitido entender con mayor detalle el 
proceso de formacion de sistemas binarios de proto-estrellas. 
Sin embargo, este es un fenomeno fisico complejo, tal que 
un estudio abordado de manera tradicional implicaria la 
imposibilidad de visualizar los detalles del proceso. Es aqui en 
donde la potencia de la informatica y las simulaciones numericas 
se vuelven una herramienta muy util, por no decir indispensable, 
a la hora de estudiar el proceso de formacion en detalle. En las 
secciones siguientes intentaremos mostrar al lector algunos de 
los aspectos fisicos mas importantes del colapso gravitacional 
de nubes progenitoras. Nos interesa destacar la inter-relacion 
de fuerzas como una manera de ilustrar la aplicacion de las leyes 
de la fisica clasica y la ocurrencia de la fragmentacion de nubes 
como el mecanismo principal para dar origen a grupos binarios 
de proto-estrellas. 



LA ETAPA INICIAL DEL COLAPSO GRAVITACIONAL 

Como se puede ver en la Figura 1, hay mucha actividad en 
las nubes progenitoras de origen tanto interna como externa. 
Por ejemplo, las nubes se calientan porque reciben radiacion 
cosmica (proveniente de regiones distantes del Universo); fuertes 
vientos estelares (provenientes posiblemente de explosiones de 
estrellas supernovas localizadas en la misma nube).Tambien es 
comun que ocurran colisiones entre nubes como consecuencia 
de fuerzas de marea provocadas por la presencia de galaxias 
vecinas o incluso por colisiones directas entre galaxias enteras, 
provocando la redistribucion del gas en la nube. 

Debido a tales eventos, en algunas regiones de la nube el 
gas se acumula en torno a pequenas perturbaciones de masa 
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que alcanzan de manera fortuita una 

densidad mayor que su medio ambiente 

circundante. La fuerza de gravedad intenta 

siempre ensamblar mas materia sobre 

estas primeras acumulaciones fortuitas de 

gas debido a su naturaleza de atraccion 

entre toda la materia inter-actuante. Sin 

embargo, la presion del gas se opone 

permanentemente a ser comprimido por 

la fuerza de gravedad. Ahora bien, se ha 

observado que las nubes progenitoras son 

en general frias: una temperatura tipica 

seria de 10 grados kelvin, es decir, unos 

263 grados Celsius bajo cero. Mientras 

que su masa puede llegar a ser tan grande 

como cientos de veces la masa de nuestro 

Sol. Asi pues, no es sorprendente que la fuerza 

de gravedad tenga muchas posibilidades de 

ganarle a la fuerza de presion en esta primera 

etapa de confrontacion entre fuerzas de tal forma que 

en varias regiones separadas de la nube se inicia el proceso de 

acumulacion de gas, ver Figura 2. 

La acumulacion inicial de gas ocurre de tal manera que una 
nube grande se ve conformada por varias subregiones de gas 
mas o menos bien definidas, de menor tamano pero de mayor 
densidad promedio que la nube mayor. Cada una de estas sub- 
nubes continua tranquilamente con el proceso de acumulacion 
de gas por gravedad durante miles de anos; aunque existe 
la posibilidad que eventos ajenos perturben el proceso de 
acumulacion. A esta clase de sub-nube se le llama nube pre- 
estelar porque ya es una fuerte candidata a iniciar el proceso de 
colapso gravitacional. 

El nacimiento de una nueva estrella tiene lugar cuando el 
colapso gravitacional de la nube pre-estelar es tan intenso que 
la temperatura supera el valor critico necesario para que tengan 
lugar los procesos de fusion nuclear propios de una estrella. 
No obstante, antes de llegar a esta etapa final del proceso de 
formacion estelar, en el transcurso del colapso de la nube pre- 
estelar, ocurren interesantes fenomenos mecanico-clasicos, 
quesepuedenvislumbraralaplicarlasleyesdelafisica 
a un modelo de nube pre-estelar simplificado, 
como hacemos a continuacion. 

LOS MODELOSDE NUBE 

Las observaciones de los astronomos 
indican que las nubes pre-estelares tienen 
nucleos mas o menos bien definidos 
que se encuentran rodeados por un halo 
envolvente de gas. El nucleo de la nube 
contiene el gas mas denso; mientras 
que el que se encuentra en las regiones 
intermedias y mas alejadas del centro de la 
nube son cada vez menos densas. A estas 
nubes se les conoce como centralmente 
condensadas. Aunque este no es el unico 
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caso observado en las nubes, porque se 
infieretambien la existencia de nubes pre- 
estelares en las que el gas se distribuye 
de manera uniforme, de tal forma que en 
cualquier region de la nube el gas tiene 
casi la misma densidad. Este tipo de nubes 
son denominadas radialmente uniformes. 
En ambos casos la densidad media de la 
nube pre-estelar es mayor que en el resto 
de la nube progenitora, razon por la cual 
el colapso gravitacional se inicia [1]. 

Como en todo fenomeno fisico, 
las condiciones iniciales determinan 
fuertemente el resultado de la evolucion 
dinamica del sistema de interes. Asi pues, 
aunque la primera parte del colapso de 
la nube se Neva a cabo de igual manera 
en cualquier modelo de nube pre-estelar, 
podemos esperar que el resultado final del mismo 
sea diferente dependiendo de si la nube pre-estelar es 
radialmente uniforme o centralmente condensada. 

Consideremos una nube pre-estelar de gas con forma esferica 
que rota en torno a un eje, con una velocidad angular constante, 
tal como se ilustra en las figuras 2 y 3. Como habiamos indicado, 
la fuerza de gravedad jala toda la materia hacia el centro de la 
nube mientras la fuerza de presion empuja el gas hacia afuera. La 
rotacion de la nube provoca la aparicion de una fuerza centripeta 
dirigida hacia el centro en cualquier piano de latitud (todo piano 
en la esfera paralelo al piano ecuatorial). De esta forma, todo 
elemento de gas sentira una aceleracion total hacia el interior 
de la nube, que esta formada por dos contribuciones: una parte 
debida a la aceleracion centripeta dirigida hacia el centro del 
piano y otra parte debida a la aceleracion gravitacional dirigida 
hacia el centro de la nube. Estas aceleraciones tienen direcciones 
co-lineales solo en el piano ecuatorial mientras que en los polos 
son perpendiculares. Por esta razon, la aceleracion total es mayor 
en el circulo ecuatorial que en los regiones cercanas a los circulos 
polares. 

Recordemos que el momento angular es una medida 

de la cantidad de rotacion de un cuerpo; ademas 

una cantidad que se conserva durante toda 

a evolucion de la nube ya que no hay torcas 

externas actuando sobre la nube. Asi, 

como consecuencia de la rotacion y por 

conservacion de momento angular, la 

nube empieza a aplanarse de los polos 

hacia el ecuador, de tal manera que se 

forma un disco aplanado en el ecuador 

que gira mas rapido que el resto del gas 

que todavia esta en la nube esferica; este 

disco se le conoce como disco de acrecion, 

porque el gas continua cayendo sobre el 

por efecto de la fuerza de gravedad. 

La anchura del disco de acrecion queda 
determinada por el equilibrio mecanico 




Figura 4. Anchura del disco de acrecion y diagrama de cuerpo libre de 
un elemento de fluido. 



concentra en la region del disco de acrecion cuya extension es de 
aproximadamente el 4 % del radio inicial. 

El mosaico de la figura 5 corresponde a un modelo de nube 
pre-estelar con densidad uniforme mientras que el mosaico 
de la figura 6 muestra los resultados del colapso para una 
nube centralmente condensada. La primera y mas importante 
diferencia entre estos mosaicos es que en la figura 5 la nube 
colapsa formando una sola concentracion de materia en el 
centro de la nube, mientras que en la figura 6 se ve que la nube 
se fragmento, de tal forma que aparecen varias concentraciones 
de materia resultantes que se encuentran orbitando unas con 
respecto a las otras. En ambos modelos se ve la formacion de 
brazos espirales alrededor de las concentraciones de materia que 
aparecen tambien como una consecuencia de la rotacion inicial 
de la nube pre-estelar. <;C6mo se explica esta diferencia en los 
resultados? 



entre las componentes verticales de las fuerzas gravitacional y de 
presion del gas, tal como mostramos en el recuadro de la figura 
4, en el cual se ve un diagrama de cuerpo libre de un elemento 
de gas. En el equilibrio hidrodinamico, la proyeccion de la fuerza 
de gravedad en la direccion vertical, debe cancelar a la fuerza de 
presion que el gas ejerce como oposicion al aplanamiento de la 
nube. 

LOS RESULTADOS DE LAS SIMULACIONES NUMERICAS 

En nuestras simulaciones numericas del colapso de nubes pre- 
estelares, hay muchas decisiones que tomar para establecer el 
estado fisico inicial del modelo de nube, a saber: hay que fijar la 
velocidad de rotacion y la temperatura; el radio y masa iniciales; 
la densidad de las nubes pueden ser diferentes y ademas, en 
el caso de nubes centralmente condensadas, la extension del 
nucleo comparado con el halo de gas tambien puede variar, 
esto es, la razon RO / RC de las figuras 3 y 4. Estas elecciones nos 
permiten investigar algun efecto en particular sobre el resultado 
final del colapso. 

Para mostrar los resultados de algunas de las simulaciones 
usamos mosaicos formados por varios paneles, como se ilustra 
en las figuras 5 y 6. Cada panel en el mosaico corresponde a un 
tiempo de evolucion del modelo, tal que el tiempo aumenta 
de izquierda a derecha y de arriba hacia abajo. Cada pequena 
imagen muestra -en una misma escala de colores para cada 
mosaico- la distribucion de densidad de una rebanada de nube 
paralela al piano ecuatorial, tal como se la vena desde arriba. Por 
ejemplo, de acuerdo con la barra de escala de densidad (que 
aparece en la parte inferior de cada mosaico) el color amarillo 
indica regiones de la nube con densidades altas; los colores 
rojo y verde indican regiones con densidades intermedias y 
por ultimo, el color azul indica regiones de baja densidad en 
la nube. Las escalas de longitud que aparecen es los ejes de 
cada panel estan normalizadas con el radio inicial de la nube y 
la escala de densidad esta normalizada con la densidad inicial 
de la nube. Se debe notar entonces, por las escalas en los ejes 
X (horizontal) e Y (el vertical), que en la figura 5 vemos el gas 
mas denso que esta contenido hasta el 1 0% del radio inicial de la 
nube, mientras que en la figura 6 vemos que el gas mas denso se 



Hemos implementado delibe-radamente una perturbacion 
de masa simetrica con respecto al origen de las coordenadas 
del piano ecuatorial de la nube, con el proposito de favorecer la 
formacion de sistemas binarios como resultado. De ahi que los 
primeros cumulos de gas que logran formarse son antipodas el 
uno del otro, de tal manera que una linea imaginaria que los une 
pasara a traves del origen de coordenadas de la nube esferica 
original. 

Cada uno de los cumulos de gas recien formados por el 
fuerte colapso gravitacional de la nube pre-estelar (mostrados 
en amarillo en las figuras 5 y 6) ejercen por lo tanto un par 
gravitacional de atraccion el uno sobre el otro. Entonces la 
velocidad de los elementos de gas empieza a alinearse con el eje 
imaginario de simetria que une los cumulos, con el efecto neto 
que este gas del cumulo pierde momento angular. Entonces 
sabemos que el gas que se acumula sobre el disco de acrecion 
tiene poco momento angular; mientras que el gas que se 
acumula en los brazos espirales tiene mayor momento angular. 
Asi es como los cumulos principals de la simulacion pierden su 
momento angular, de tal forma que la fuerza de gravedad que 
se ejercen entre ellos, las hace acercarse cada vez mas hasta 
que finalmente se fusionan, como se puede ver en los primeros 
paneles de la figura 5. 



O10H 
-0.00 I 



0.10 1 

ol:-i 
-o.oo 




o.io ■ 
-o.oo I 



oo? ■ ■ OtoH ■ ocoH 

"0-ioH I oiuH ■ oi[»H 

o.io o.tD oou g.[» u.Lo g.ui uju^ o.oo u.oi o.io u.L» o.o> o.cio og> ».u 

10 

o:i-> 

o:ki 

0.05 ■ ■ -0.0? ■ ■ -0.05 1 

OKIE ■ OILlR ■ IL>I 

o.io uo5 oou u.cs g.lo u.ui uoi o.oo ooi 0.10 u.Lo o.o> am aa\ ftll 






CHI-,9 
CHKlH 
-0.05 H 

-oioH 



*z 



-0.05 ■ ■ OtoH ■ -0.05 

oio^H^^^H^H -o.io^^^^HH -o.iol 

-0.lO-0.G5-CS.UGG.CKi U. LO -G.1CI -GOi -0.00 U.GJ O.IO -G.lO 




(j la r 

005 1 



0.10| 
0.05 
-G.OCll 




I 0,05 ■ 
■ OLKlH 
omH I OL!>.|| ■ 0(11 ■ 

■■■ l:: I - -M -O.iol 

0.10 0.05 O.OO 0.05 0.L0 0.10 0.05 0.00 0.05 O.IO 0. LO 0.05 0.00 0.05 It 



Figura 5. Mosaico de una simulacion con una sola proto-estrella central 
como resultado. 
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Figura 6. Mosaico de una simulacion con varias proto-estrellas como 
resultado de la ocurrencia de fragmentation. 

Ahora bien, la causa por la que no se produce fusion de 
cumulos en la simulacion de la figura 6, es tal vez la existencia 
de pequenas variaciones en las posiciones y velocidades de 
las particulas que forman el gas, que tienen su origen en la 
aleatoriedad de la distribucion de particulas que representan a la 
nubepre-estelar inicial. 

iQue determina que ocurra o no fragmentacion de la nube 
pre-estelar durante su colapso? 



LA FRAGMENTACION. 

Para que una teoria de formacion estelar tenga posibilidad de 
ser exitosa, debe explicar las observaciones astronomicas que 
muestran la presencia de miles de grupos binarios y multiples de 
estrellas formadas en el seno de una sola nube progenitora. Se 
necesita de algun mecanismo fisico que provoque la ocurrencia 
de la fragmentacion de la nube pre-estelar, algo como se ilustra 
esquematicamente en la figura 7. 




Figura 7. Diagrama esquematico que representa el proceso 
de fragmentacion teorico requerido por las observaciones 
astronomicas. 

Fue Fred Hoyle, astrofisico ingles autor de la Teoria del Estado 
Estacionario, el primer cientifico que propuso una explicacion 
para la fragmentacion durante el colapso gravitacional de 
la nube [2] . Hoyle argumenta que en la medida en que la 
nube colapsa, los grumos de materia (necesarios para iniciar 



nuevas aglomeraciones de materia que pueden veneer la 
presion termica y promuevan a su vez nuevamente el colapso 
gravitacional) pueden ser cada vez menores en la medida en 
que la temperatura de la nube no aumente por efecto del propio 
colapso. Sin embargo, sabemos que en la medida en que la nube 
colapsa, su tamano se reduce; y que por esta reduccion del radio, 
la energia gravitacional se convierte en calor. 

iComo se puede mantener la nube con una temperatura 
aproximadamente constante durante el colapso? Cuando el 
colapso no es todavia muy intenso, el gas se puede liberar de este 
exceso de energia cinetica mediante emision de radiacion de las 
moleculas excitadas de hidrogeno. La radiacion emitida debe 
ser capaz de difundirse hasta la superficie exterior de la nube y 
escapar finalmente para que este mecanismo de enfriamiento sea 
efectivo y permita la continuacion del colapso isotermico, el cual 
podria tener como consecuencia la tan deseada fragmentacion. 
No obstante, en la medida en la que la densidad de la nube 
aumenta por efecto del colapso, la difusion de la radiacion a 
traves del medio molecular se hace cada vez mas dificil. Entonces, 
como era de esperarse, el aumento de la opacidad de la nube 
hace que la temperatura del gas empiece a aumentar. En otras 
palabras, el aumento de temperatura hace que el gas aumente 
su resistencia al colapso dotandolo de mayor presion termica. 
Llega un momento en que incluso el colapso se hace mas lento 
(o se detiene por completo) y la fragmentacion es mucho menos 
probable que ocurra, ya que no habra varias regiones separadas 
en la nube en las que puedan crecer perturbaciones de masa 
simultaneamente. 

En el trascurso de nuestras investigaciones, hemos observado 
que no solo son los efectos termodinamicos de la nube los 
que determinan la ocurrencia de fragmentacion. Tambien hay 
elementos de origen puramente dinamicos que tienen influencia 
sobre la fragmentacion. 

De hecho, en las referencias 4, 5 y 6 hemos reportado varias 
maneras en que la fragmentacion de las nubes pre-estelares 
puede ocurrir. Una manera es que la barra central se fragmente; 
otra seria que los brazos espirales se rompan y se separen del 
cumulo central debido a que sus velocidades de rotacion son 
diferentes; tambien pueden ocurrir ambas posibilidades en un 
mismo modelo, como se ilustra en la Figura 8. En esta figura se 
muestra que se produce la fragmentacion de la barra central 
dando lugar a la formacion de dos cumulos proto-estelares 
interiores. Cabe senalar que los dos fragmentos exteriores 
resultantes del rompimiento de los brazos espirales estan 
ya presentes en el momento en que la fragmentacion de la 
barra central se produce. En este momento, el resultado de la 
simulacion son cuatro fragmentos proto-estelares. Poco despues, 
la ocurrencia de fusion entre dos fragmentos hace que solo 
quedan dos fragmentos que entran finalmente en orbita uno con 
respecto al otro hasta alcanzar el equilibrio hidrodinamico, como 
se ilustra en la Figura 9. 

En el articulo de investigacion de la referencia 3, encontramos 
que a medida que la region envolvente del nucleo de la nube 
aumenta en extension (esto es, la razon R°/R c aumenta), la 
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Figura 8. Acercamiento 
a la region central de 
un modelo con multiple 
fragmentacion. 

posibilidad de un modelo para 

producir fragmentos proto-estelares 

disminuye. Para explicar por que 

ocurre esto, conviene enfatizar 

que todas las nubes con diferentes 

razones R° / Re fueron dotadas 

inicialmente de la misma energia 

rotacional. Pero en el caso de las 

nubes centralmentecondensadas,esta 

energia rotacional debe compartirse 

entre el nucleo y el gas envolvente. 

Entonces, mientras mas grande sea el 

tamano de la capa envolvente mayor energia 

debe tomar de la nube para mantenerse rotando 

con la misma velocidad angular que el nucleo; mientras 

que al nucleo le queda menos energia rotacional disponible para 

presentar oposicion al colapso gravitacional; entonces, mientras 

mayor es el envolvente se induce un colapso mas fuerte en la 

region central que provoca la destruccion de cualquier estructura 

que se haya formado durante las fases iniciales del colapso. Asi, 

dotar al nucleo de suficiente energia de rotacion es crucial para 

que ocurra la fragmentacion. 



CONCLUSIONES. 

En el curso del colapso, la densidad de la nube pre-estelar 
aumenta en un rango que va desde 10~ 18 gr/cm 3 hasta 10~ 8 gr/ 
cm3, en la cual ya podemos identificar a los cumulos de gas 
mas densas de la simulacion como proto-estrellas. Se sabe que 
muchas de las caracteristicas dinamicas de una proto-estrella 
seran heredadas por la estrella real que se formaria de seguir el 
proceso de colapso, hasta densidades del orden de 10~ 1 gr/cm 3 . 

Nosotros nemos logrado seguir la evolucion del colapso hasta 
estas densidades intermedias. En este articulo de divulgacion 
hemos presentado brevemente algunos aspectos importantes 
del colapso de nubes pre-estelares. Una de nuestras conclusiones 
es que el resultado final del colapso de nubes, y en particular 
la ocurrencia de la fragmentacion, depende fuertemente de 
aspectos fisicos tanto termodinamicos como dinamicos. 

Finalmente, deseamos compartir con el lector 
nuestra conviccion de que el paradigma 
tradicional de la investigacion en fisica (teoria 
y experimento) se ha visto enriquecido 
con el desarrollo de la computacion. 
Hoy dia los investigadores estamos en 
posicion de usar la enorme capacidad 
de calculo de los computadores 
actuales para simular procesos 
fisicos muy complejos como el de 
formacion estelar. Pensemos que los 
resultados de una simulacion que 
se conllevan con las observaciones 
astronomicas, es un claro indicativo 
de que el modelo elegido se 




encuentra en el buen camino, pero aun 
cuando estos se desvian claramente 
de los datos observacionales siguen 
siendo utiles poniendo en evidencia 
que el modelo es erroneo o por lo 
menos incompleto, pudiendo haber 
variables importantes que aun no 
han sido tenidas en cuenta. 
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IDENTIFICACION DE IMAGENES Y VOCES MEDIANTE 

REDES DE HOPFIELD 

UNA APLICACION DE LAS REDES NEURONALES RECURRENTES 



ARNULFO CASTELLANOS MORENO, FABIOLA SUSETH LOPEZ AGUIRRE 




^~ n la Universidad de Sonora podria desarrollarse un sistema de reconocimiento de voces utilizando 
^™ una red neuronal como la que sera explicada al final de este articulo. Lo que se requiere es el sistema 
^— electronico para reconocer las frecuencias mas importantes de cada voz, medir la intensidad del armonico 
correspondiente y guardarla en un arreglo numerico del tipo de los vectores que se aprenden en algebra 
lineal. A falta de un dispositivo similar, ejemplificamos aqui con el reconocimiento de letras, sin embargo, el 
procedimiento seria fundamentalmente el mismo. 
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jQUE SON LAS REDES NEURONALES 
RECURRENTES? 

Las redes neuronales recurrentes son sistemas ideados 
como abstracciones de las estructuras neurobiologicas 
(cerebros) encontradas en la naturaleza y tienen la 
caracteristica de ser sistemas desordenados capaces 
de guardar informacion. La forma en que desarrollan 
su trabajo es esencialmente distinta de la utilizada por 
las computadoras convencionales. Los procesadores 
microscopicos del cerebro (neuronas) operan en paralelo 
y presentan cualitativamente mas ruido que los elementos 
que forman las computadoras. No ejecutan un programa 
fijo con base en un conjunto previamente especificado 
de datos, sino que comunican senales a traves de 
retransmisores que llamamos sinapsis, que Megan a 
centros de conjuncion llamados cuerpos de las neuronas 
y desde los cuales surgen senales electricas a traves de 
canales conocidos con el nombre de axones. 

La importancia de cada sinapsis en el proceso de 
retransmision se actualiza continuamente y lo mismo 
ocurreconalgunaspropiedadesintrinsecasde las neuronas, 
proporcionando un sistema de autoprogramacion y 
adaptacion que sustituye a la programacion externa de los 
sistemas de computo comunes. Existe asi una dinamica de 
las sinapsis y de las neuronas en el cual los programas y los 
datos cambian todo el tiempo. 

Los cerebros no tienen una arquitectura o mapa de 
conexiones especifico, por ejemplo, presentan estructuras 
relativamente ordenadas en centros receptores de 
informacion, pero en las regiones procesadoras de 
alto nivel como el neocortex, son estructuras amorfas, 
practicamente aleatorias. 

Operando en paralelo, los cerebros dedican grandes 
cantidades de neuronas a una actividad determinada, 
lo cual da por resultado muchas actividades posibles, 
diferentes y complejas. Pueden adaptarse a circunstancias 
cambiantes, reconocer objetos deformados, borrosos, o 
parcialmente visibles. En smtesis, las redes neuronales son 
esencialmente diferentes de los sistemas de computo de 
la vida cotidiana. 

Los distintosmodelosde redes neuronales hanplanteado 
importantes e interesantes retos a los fisicos teoricos, pues 
exhiben conductas muy ricas y nada triviales. Se trata de 
sistemas sumamente complejos de analizar por el gran 
numero de elementos interactuantes, por el caracter no 
lineal de los elementos individuals que operan entre ellos, 



porque las interacciones no necesariamente son identicas 
ni simetricas y porque la magnitud de las interacciones, 
asi como los umbrales internos de cada neurona, pueden 
cambiaren el tiempo. 

En los intentos por sistematizar la comprension del 
cerebro pueden distinguirse dos caminos: i) el primero es 
el seguido por los neurofisiologos, quienes generalmente 
confrontan una abundancia de datos tan vasta que 
dificulta el examen sistematico dentro de un marco de 
trabajo experimental. En este enfoque, los datos estan 
perturbados por una cantidad muy grande de influencias 
que no pueden ser eliminadas facilmente. ii) El segundo 
es el que han abierto los constructores de modelos sobre 
la base de un conjunto tambien grande de suposiciones 
acerca del sistema que se busca explicar, con la esperanza 
de que el enfoque ayude al surgimiento de nuevas 
hipotesis, susceptibles de ser verificadas con la diversidad 
de datos disponibles. 

Dentro del segundo de los caminos hay dos opciones 
para los creadores de modelos: 1°) la que consiste en tratar 
de modelar al sistema real tanto como sea posible, pero 
con el dilema de que se plantean tantos parametros que 
en realidad no se alcanza un conocimiento profundo; 2°) 
la opcion de descartar, a priori, todos aquellos parametros 
que a primera vista parecen no ser esenciales, a fin de 
simplificar el analisis matematico. 

Los primeros buscan tanta fidelidad al sistema real 
que terminan siendo una copia de mala calidad en la 
que la comprension del fenomeno se desvanece. Los 
segundos no constituyen representaciones realistas 
del cerebro, sino que su inspiracion neuronal puede 
contribuir a comprender algunas de las propiedades que 
los caracterizan en el procesamiento de informacion. Entre 
otras, de aqui surgen las redes neuronales recurrentes. 

Tomando como opcion a la segunda de las mencionadas 
arriba, se ha definido lo que se llama una neurona formal 
[1], y sobre la base de esta, a la red neuronal en su forma 
mas general. El primer paso consiste en definir como 
sigue los estados de una neurona: a) el estado de lanzar 
una serial, y b) el estado de inactividad. Al primero se le 
asigna el numero +1 y al segundo el -1. El segundo paso 
es establecer la conexion entre las neuronas mediante una 
matriz {J.j} en la que se especifica cuanto vale la conexion 
de la i-esima con la j-esima neurona. El tercer paso es 
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establecer una relacion funcional en la que la j-esima 
neurona recibe la informacion de todas las demas y en la 
que dependiendo de su umbral, produce, o no produce, 
una serial de intensidad F. . En terminos matematicos la 
forma es 




Para que la red neuronal pueda operar, es conveniente 
que las funciones Fi sean estrictamente crecientes y 
acotadas, del genero de las llamadas funciones sigmoides, 
como es el caso de la tangente hiperbolica y de la funcion 
signo que se presentara mas adelante. Aqui es donde se 
introduce la no linealidad mencionada anteriormente y 
tambien el umbral de cada neurona, que es el limite de 
recepciones que puede soportar antes de mandar una 
serial. 

Dependiendo del tipo de red, los pesos son variables 
que seran modificadas durante su evolucion, o por el 
contrario, seran magnitudes fijas, como en el modelo de 
Hopfield que se discute en la siguiente seccion. 

Las neuronas tienen entradas y salidas, tal que estas 
ultimas son calculadas con base en las primeras. Asi surge 
el concepto de regla de aprendizaje, que consiste en la 
sistematizacion de la asociacion entre una pregunta y una 
respuesta. En general, la pregunta r es un vector de mayor 
dimension que la respuesta R. Un ejemplo tipico es el 
reconocimiento de voces o sonidos, en los cuales el vector 
de entrada r es un espectro de frecuencias, que debe ser 
transformado en una lista de palabras R. 

UN EJEMPLO DE RED RECURRENTE ES LA DE HOPFIELD. 

Un ejemplo simple de como se hace esto es la siguiente 
red de Hopfield simplificada [2]: 

Consideremos una red de n neuronas cuyo estado esta 
dado por un vector r=(x 1 x 2 „..., x n ), que nos indica si la 
i-esima neurona esta enviando serial con un+1,ocon un-1 
si esta apagada. La operacion de la red la podemos imaginar 
como si fueran las luces de un arbol navideno en los que 
cada foquito se prende y se apaga. Si en determinado 
momento tomamos una fotografia con una camara de alta 
resolucion y alta velocidad, lo que habremos logrado es 
captar uno de esos vectores r. Ademas se necesita el peso 
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de la conexion w.. de la neurona i con la neurona j. Tambien 

u J 

se necesitan las memorias a guardar, que podrian ser, por 
ejemplo, el catalogo de las frecuencias de las voces con las 
cuales sera comparada la voz que se desea identificar. Este 
catalogo es un conjunto de p memorias que denotaremos 
como:{0u} 

Los pesos de la conexion estaran dados por la siguiente 
matriz 



™ij = 



! 



si i = j 



que puede ser organizada en un arreglo rectangular en 
la forma siguiente: 
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La capacidad de las neuronas para soportar la recepcion 
de senales, sin transmitir ellas, se introduce mediante los 
umbrales de las neuronas organizados en el vector que 
sigue:e = (6 1 ,e 2 ,...,e n ) 

La decision de emitir, o no emitir, se especifica con 
la regla de actualizacion del vector de estados de la red 
neuronal. Esta es 

r(K+l)=sign[wirck)-9] 

Para saber cuando el estado alcanzado por la red 
puede ser considerado como satisfactorio, se especifica 
un regimen de tolerancia tal que, una vez cumplida la 
siguiente desigualdad 

\xi(k + l) -Xj(k)\ <s 

la red deja de operar. 

Enseguida se imprime el estado final: r fina| 
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Figura 1: una letra formada sirve como punto de partida 
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Figura 2: uno de los estados de la red cuando avanza hacia 

el reconocimiento. 
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Figura 3: un estado de la red cercano al final, en que la letra 



Figura 4: Resultado final cuando la letra es identificada. 



ya esta muy cerca de ser reconocida. 

UN EJEMPLO MEDIANTE RECONOCIMIENTO DE 
LETRAS. 

Los patrones, o memorias {0u} que se guardan en la 
red neuronal pueden ser tornados de la informacion 
codificada de un conjunto de imagenes, o de frecuencias 
de voces, etcetera. Aqui desarrollaremos un ejemplo de 
reconocimiento de letras. 

Se escoge una letra del alfabeto y se dibuja sobre una 
cuadricula. Enseguida se localizan los cuadros que tienen 
negro debido a la letra para codificarlos mediante 1-s, a su 
vez, los cuadros que tienen bianco se les asignan 0-s. Este 
es uno de los patrones o memorias a almacenar. Enseguida 
se toman mas letras del alfabeto y se hace lo mismo hasta 
completar los p patrones. Los estados s de las neuronas 
se modifican mediante la regla algebraica 2s-1 para tener 
+1-s y -1-s. 

El estado inicial se toma especificando cualquier letra 
deformada para que se vea borrosa, o parecida a otra, por 
ejemplo la E y la L, o bien la O y la U. 

Se opera la red con el algoritmo arriba mencionado y se 
obtienen resultados similares a los del siguiente ejemplo: 

Partimos de una letra X deformada como se indica en la 
figura siguiente. A la izquierda tenemos la matriz de 1-s y 
0-s que sirven para alimentar a la red neuronal, ademas de 
otras letras que se cargan en su memoria. A la derecha se 
observa la letra deformada. 



CONCLUSIONES 

La red de Hopfield que hemos descrito es un sistema 
sencillo que se inspira en las estructuras de las redes de 
neuronas biologicas y tiene aplicaciones practicas en 
el reconocimiento de imagenes, voces y senales cuya 
estructura es conocida. Puede tener aplicaciones en 
criminologia, en medicina para ubicar enfermedades 
potenciales a partir de cuadros de diagnostico, etcetera. 
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La contaminacion de suelos es un problema recurrente en muchas regiones del mundo. Cuando los 
residuos organicos 6 inorganicos, producidos por ciertas actividades humanas, son derramados en el 
suelo y afectan la calidad de vida de los seres vivos, se consideran contaminantes segun la normatividad 
de proteccion ambiental de nuestro pais. En este trabajo se describe una metodologia teorico experimental 
para entender los fenomenos que ocurren durante el proceso de eliminacion de contaminantes en suelos 
aplicando campos electricos, a fin de optimizar el rendimiento de dicho proceso. Esta tecnologia de punta 
actualmente es estudiada en algunos paises desarrollados. 
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EL PROBLEMA GLOBAL DE LA CONTAMINACION DE 
SUELOS 

Historicamente, los suelos de nuestro planeta han sido 
contaminados como resultado de diversas actividades humanas, 
tales como la agricultura (1 ) el tratamiento de la madera, la extraccion 
y procesamiento de metales y los conflictos belicos, entre otras (2). 
En general, los contaminantes de suelos se pueden clasificar en dos 
grandes grupos: organicos e inorganicos. Dentro de los primeros 
se tienen los aceites y gasolinas provenientes de la industria del 
petroleo, los fertilizantes, fungicidas, plasticos, etc. Dentro de los 
segundos se clasifican una variedad de compuestos quimicos que 
contienen metales y otros elementos, tales como el arsenico, cobre, 
plomo, cadmio, mercurio, etc. A estos ultimos se les denomina 
comunmente metales pesados. Recientemente, la contaminacion 
por metales pesados en sitios de minas y fabricas abandonadas ha 
causado serios problemas al medio ambiente y a la salud humana. 
Como ejemplossetiene la contaminacion por fluorita enCorea(3) y 
por plomo en Mexico (4). 

Adicionalmente al dano ecologico que los suelos contaminados 
ocasionan al habitat de un determinado lugar, la contaminacion 
de los suelos cercanos a las ciudades representa perdidas de 
oportunidades para la instalacion de nuevos centros de trabajo, y por 
lo tanto para el desarrollo economico. Por lo tanto, el uso de estos 
sitios como receptores para confinar residuos peligrosos no es una 
opcion. Limpiarlos para darles un mejor uso es una mejor alternativa 
la cual es congruente con la doble responsabilidad de cuidar el 
medio ambiente y propiciar el desarrollo social y economico de las 
poblaciones. 

LIMPIEZA DE SUELOS POR MEDIO DE ELECTRICIDAD 

Existen varios metodos para limpiar suelos contaminados. Este 
articulo se centra en el proceso de eliminacion electrocinetica, 
cuya tecnologia novedosa es foco de atencion actual en los paises 
desarrollados, tales como Estados Unidos y Holanda (5). El principio 
de funcionamiento de este proceso se muestra en la Figura 1 y 
consiste en aplicar un campo electrico al suelo contaminado a 
traves de una serie de electrodos para provocar que los compuestos 
metalicos disueltos en el agua del suelo se muevan por efecto del 
campo electrico hacia alguno de los electrodos, dependiendo 
de la carga electrica del compuesto metalico en cuestion. Asi, los 
compuestos metalicos con carga positiva se desplazaran hacia el 
catodo, y aquellos con carga negativa lo haran hacia el anodo. Al 
cabo de un tiempo determinado, los contaminantes se concentraran 
en areas localizadas del suelo, de donde podran ser extraidos para su 
disposicion final (6, 7). 

El proceso de eliminacion electrocinetica comprende 
fundamentalmente tres fenomenos de transporte que ocurren 
simultaneamente. El primero es causado por el campo electrico 
que hace que los iones metalicos (atomos o grupos de atomos con 
carga electrica) se muevan a traves del liquido contenido en los 
poros del suelo hacia el electrodo con carga opuesta a la del ion. 
A este fenomeno se le denomina electromigracion. El segundo 
fenomeno ocurre muy cerca de la superficie de los poros del suelo. 
La superficie de las particulas que forman el suelo tipicamente 
posee carga electrica. El liquido cercano a la superficie adquiere 
una carga electrica de signo contrario a la carga de la superficie de 
los poros. De nueva cuenta, este liquido es afectado por el campo 
electrico y se mueve hacia alguno de los electrodos. Arrastrando 
consigo el resto del liquido contenido en los poros por efectos 



viscosos. El movimiento del liquido se Neva los contaminantes que se 
encuentran dentro de el. A este fenomeno se le llama electroosmosis. 
Un tercer fenomeno se presenta cuando en una region del suelo los 
contaminantes se encuentran con una concentracion mayor a la 
de otra region del mismo suelo. En este caso, los contaminantes se 
mueven espontaneamente de la region de mayor concentracion a 
la region de menor concentracion. A este fenomeno se le conoce 
como difusion ordinaria y no requiere de la presencia de un campo 
electrico. 

Simultaneamente a los fenomenos de electromigracion, 
electroosmosis y difusion ordinaria descritos anteriormente, ocurren 
reacciones quimicas entre los diferentes compuestos disueltos en el 
liquido del poro, entre estos y las particulas de suelo. Adicionalmente, 
en los electrodos negativos (catodos) y positivos (anodos) se 
presenta la descomposicion electrolitica del agua. Esta consiste en 
la desaparicion de las moleculas de agua por accion de la corriente 
electrica para producir, respectivamente, hidrogeno y oxigeno los 
cuales se desprenden hacia la atmosfera en forma de gas. En los 
anodos tambien se generan electrones, los cuales son impulsados a 
traves de la fuente de poder hacia los catodos. 
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Figura 1. Principio de funcionamiento del proceso de eliminacion 
electrocinetica 

La eliminacion electrocinetica es un proceso de bajo costo y 
eficiente (8) para descontaminar suelos con metales pesados. Sin 
embargo, antes de llevar a cabo su aplicacion en gran escala es 
necesario realizar experimentos en el laboratorio con el fin estudiar 
sus caracteristicas y optimizar su funcionamiento bajo condiciones 
controladas, en funcion del contaminante y tipo de suelo a tratar. 
Esta es una etapa decisiva en el desarrollo del proceso, ya que en 
cuestion de limpieza de suelos, como en muchos otros campos de 
la ingenieria, no existen panaceas, es decir formulas magicas de uso 
universal. 

La forma mas comun de estudiar este proceso es mediante el uso 
de prototipos llamados celdas experimentales a nivel laboratorio, 
como la mostrada en la Figura 2. La celda consiste de un tubo de 
acrilico o vidrio dentro del cual se coloca el suelo o medio poroso con 
el contaminante. El tubo se cierra por ambos extremos y en estos se 
colocan los electrodos. Posteriormente se satura el suelo dentro del 
tubo con agua. Mediante experimentos controlados, el investigador 
aplica un voltaje conocido entre los extremos de la celda y conforme 
transcurre el tiempo se monitorea la concentracion del contaminante 
y el pH del suelo a lo largo del tubo. A partir de esta informacion es 
posible determinar el grado de eficiencia del proceso. 

EPISTEMUS 43 



Campo electrico 







Fuente de poder 














n 




J_ 


+ 

o 

c 

< 


^r^Muestra de suelo jgr^ 


■g 



Figura 2. Esquema del sistema experimental del proceso de 
eliminacion electrocinetica en una dimension. 



A la fecha, la mayoria de los trabajos sobre eliminacion 
electrocinetica se nan limitado a estudios en una sola dimension, 
como el mostrado en la Figura 2. Esto implica una serie de 
desventajas. Por ejemplo, solamente se pueden cuantificar los 
fenomenos que ocurren entre dos electrodos; adicionalmente, 
el sistema no representa la geometria de un suelo real. Por lo 
tanto, realizar experimentos en dos dimensiones con mas de 
dos electrodos representaria un avance, ya que el sistema se 
asemejaria en mayor grado a la realidad. 

METODOLOGIA DE ESTUDIO Y EVALUACION 

La Figura 3 representa la metodologia desarrollada por lo 
autores, la cual consta de dos partes: una parte teorica que 
consiste de un modelo computacional para representar los 
fenomenos de electromigracion, electroosmosis y difusion 
ordinaria descritos con anterioridad, y una parte experimental. 
El modelo matematico es una adaptacion del codigo original 
desarrollado por Jacobs y Probstein en el Massachuchets 
Institute of Technology (6, 9). La parte experimental consiste de 
una celda a escala laboratorio en forma de triangulo equilatero 
con las esquinas redondeadas, la cual se muestra en la Figura 4. 
La forma triangular es con el fin de asegurar la geometria minima 
para un estudio en dos dimensiones; es decir, en forma de un 
piano horizontal. Acorde a las condiciones del Estado de Sonora, 
como ejemplo se considera la eliminacion de cobre como 
contaminante principal. Para simular el suelo se usan particulas 
de silica gel de tamano uniforme (846 micras). 

Sobre este sistema se inserta un conjunto de electrodos 
pasivos para registrar los valores de voltaje una vez que se aplica 
una diferencia de potencial entre el anodo y los catodos. Los 
electrodos pasivos se localizan en cada uno de los puntos de 
muestreo mostrados en la Figura 4, en los cuales se mide el pH 
y se recolectan muestras de solucion acuosa para determinar 
la concentracion de cobre. Mediante un programa de computo 
desarrollado por los autores, los datos son capturados 
automaticamente y almacenados en una computadora. 
Posteriormente son procesados, visualizados graficamente 
mediante un software especializado (10) para su posterior 
interpretacion. Paralelamente se llevan a cabo simulaciones con el 
modelo computacional (6, 9) y se comparan sus predicciones con 
los datos experimentales obtenidos en la celda bidimensional. 
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Figura 3. Metodologia general para el estudio del proceso de eliminacion 
electrocinetica 



Anodo 




Figura 4. Esquema de la distribucion de los puntos de muestro para pH, 
solucion de Cu y electrodos pasivos en la celda experimental. 

RESULTADOS OBTENIDOS 

Los datos que se muestran a continuacion representan el promedio 
de varios experimentos en los cuales se colocaron particulas de 
silica gel (Si02) en la celda experimental y se empacaron con una 
maquina compactadora electrica. Este dispositivo es indispensable 
para asegurar la reproducibilidad de los datos experimentales, y fue 
disenada ex profeso para las dimensiones y geometria de la celda. 
Una vez compactado, el medio poroso se saturo con una solucion 
acuosa de sulfato cuprico para asegurar una concentracion inicial de 
cobre de 1 00 mg/L y simular el suelo contaminado. Para el inicio de la 
operacion se aplico un voltaje de 20 V entre electrodos activos (anodo 
y dos catodos) y conforme transcurrio el tiempo se registraron el pH, 
concentracion de Cu y voltaje en diversos puntos del medio poroso. 

Por brevedad, se presentan los resultados al cabo de un tiempo 
de 96 horas. La Figura 5 muestra los contornos de voltaje medidos 
experimentalmente en la celda bidimensional (5a) y los calculados 
por el modelo computacional (5b). Considerando la complejidad 
del sistema bajo estudio, la comparacion entre ambas figura 
indica que el modelo computacional reproduce razonablemente 
el comportamiento experimental, y puede ser utilizado como 
herramienta para posteriores estudios de evaluacion y optimizacion 
del proceso. 
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la) (b) 

Figura 5. Potencial electrico en el arreglo triangular al cabo de 96 horas. 

(a) Valores experimentales; (b) Valores calculados con el modelo 
computacional. 



La Figura 6 muestra los contornos de concentracion de cobre total 
en la fase acuosa del medio poroso. En este caso el modelo predice 
que el movimiento del contaminante hacia los electrodos ocurre a 
una rapidez ligeramente mayor a lo observado en los experimentos. 




N 



Figura 6. Concentracion de cobre en el arreglo triangular al cabo de 96 
horas; (a) Valores experimentales; (b) Valores calculados con el modelo 
computacional. 

Similarmente, la Figura 7 muestra los valores de pH. Se puede 
observar que alrededor de los catodos los valores de pH se 
incrementan pronunciadamente. El area de pH calculado por 
el modelo matematico es ligeramente mayor en extension y en 
magnitud que los valores experimentales. 





m 



Figura 7 Valores de pH en el arreglo triangular al cabo de 96 horas; (a) Valores 
experimentales; (b) Valores calculados con el modelo computacional. 

Las discrepancias entre los valores experimentales y las 
predicciones teoricas mostrados en las Figuras 5 a 7 motivan al 
perfeccionamiento continuo del modelo. De acuerdo a la Figura 
6a, alterminode 96 horas alrededor del 32% del area total ocupada 
por el medio poroso ha disminuido su concentracion de cobre 
hasta valores inferiores a 50 mg/L, que es el limite de tolerancia 
en suelos. Potencialmente es posible mejorar el rendimiento del 
proceso incrementando el potencial electrico aplicado entre 
electrodos e incorporando agua de lavado para mantener un 
flujo constante de liquido y extraer el contaminante del medio 
poroso. Experimentos adicionales pueden realizarse a este 
respecto en la celda bidimensional para verificar estas hipotesis. 
Para campos de suelos extensos en los cuales la experimentacion 
no es factible o muy costosa, una opcion alternativa consiste en 
realizar predicciones con el modelo computacional con el fin 
de evaluar anticipadamente diferentes opciones de operacion, 
a partir de las cuales sea posible establecer una estrategia de 
limpieza del suelo en cuestion. 

CONCLUSIONES 

Se presenta una metodologia desarrollada en el departamento 
de Ingenieria Quimica y Metalurgia de la Universidad de Sonora 
para el estudio y evaluacion de la eliminacion electrocinetica, 



el cual es un proceso novedoso para la limpieza de suelos 
contaminados por medio de electricidad. La metodologia 
consta de un modelo computacional y un sistema experimental 
bidimensional, los cuales son utilizados de manera simultanea 
para esclarecer y cuantificar los fenomenos relevantes del 
proceso, asi como para evaluar su funcionamiento bajo distintas 
condiciones de operacion. Con el debido perfeccionamiento de 
ambos sistemas, sera posible en el futuro planear y optimizar 
operaciones de limpieza de suelos reduciendo costos de 
operacion. 
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UNA APLICACION DELTEOREMA DE BAYES EN PRUEBAS DE 

DIAGNOSTICO MEDICO. 



ALEJANDRINA BAUTISTA JACOBO 



Se presenta una aplicacion del Teorema de Bayes en el area de medicina, particularmente en el campo de 
las pruebas de diagnostico. Es de gran importancia para el medico determinar si un paciente realmente 
presenta una enfermedad dado que una prueba de diagnostico le esta dando un valor positive o en 
su defecto, determinar que el paciente no presenta ninguna enfermedad alarmante debido a que se esta 
basando en un resultado negativo de una prueba. En las siguientes lineas se desarrollaran las formulas para 
determinar este grado de certeza utilizando para ello el Teorema de Bayes. 
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LA APORTACION DE THOMAS 
BAYES 

Thomas Bayes (1702-1761) fue un 
matematico britanico que aporto 
mucho al campo de la probabilidad. 
Particularmente, el teorema que 
Neva su nombre y que refiere 
a la probabilidad de un evento 
condicionado por la ocurrencia 
de otro, fue utilizado en el campo 
de la medicina, en especial para 
las pruebas de diagnostico. Una 
prueba diagnostica permite detectar 
la presencia o la ausencia de una 

enfermedad [1]. Para ilustrar esta aplicacion, supongase que 
se tiene un estudio de casos y controles, donde el evento E son 
los enfermos (casos) y S los sanos (controles) representados 
estos dos eventos mutuamente excluyentes en una tabla de 
contingencia de 2x2. De igual forma sea R+ el evento de dar un 
diagnostico positivo a la enfermedad en una prueba y R- el de dar 
un resultado negativo.Ver tabla 1 




De modo que la probabilidad 



Resultado 
de ta prueba 


Erfemno 

E 


Sano 
S 




R+ 


P{E-R+) 


P(S-R+) 


P(R + ) 


R~ 






P.R- 




P(E) 


P(S) 





Tabla 1 Presenta los posibles resultados en 
una prueba de diagnostico y las probabilidades 
marginales 



En la tabla se muestran las probabilidades marginales P(E) 
que representa la probabilidad de estar enfermo, P(S) es la 
probabilidad de estar sano, P(R+) es la probabilidad de dar un 
resultado positivo en la prueba y P(R-) la probabilidad de dar 
un resultado negativo. Estas probabilidades regularmente se 
conocen, cosa que no sucede con las probabilidades conjuntas 
P(En R+) y P(S HR+). La primera representa la probabilidad de 
estar enfermo y dar un resultado positivo en la prueba y la 
segunda indica la probabilidad de estar sano y dar un resultado 
positivo en la prueba. 

Particularmente el medico esta interesado en conocer la 
probabilidad de que el paciente presente la enfermedad 
dado que la prueba arroja un valor positivo. Esta probabilidad 
condicional se puede expresar como P(E/R+). 

Sabemos que la probabilidad condicional de un evento A dado 
que ha ocurrido B se define [2] como: 



De igual forma se cumple que 

P(Ac\B) 



1 P(A/B) = 



P(B) 






P(B/A) = 


P(AnB) 

P(A) 


Duscada se puec 


le expresar como: 


P(E/R+) = 


P{Er^R+) 
P(R+) 



De la tabla 1 se puede observar que: 

P(EnR+) 



P(EJR+): 



P(EnR+)+P(SnR+) 

Si sustituimos esta expresion en la ecuacion (1) tenemos que: 

P(Er>R+) = P(R+)*P(EfR+) 

De la ecuacion (1 ) podemos despejar el termino del numerador: 

3 P(EnR+) = P(E)*P{R+iE) 



Sabemos que P(E n R+)=P(R+ n E) [3] y podemos ver que 
tambien se cumple: 

Analogamente tenemos que: 

4 P(Sr^R+) = P(S)*P(R + /S) 

Si sustituimos las ecuaciones (3) y (4) por su termino equivalente 
en la ecuacion (2), tenemos que: 



5 P(E/R+): 



P(E)*P(R + /E) 



P(E)*P(R+fE) + P(S)*P(R+fS) 



La ecuacion resultante de estas sustituciones es llamada el 
Teorema de Bayes. 

VALOR PREDICTIVO POSITIVO (VPP) 

Ahora bien, en medicina, la probabilidad de estar enfermo 
dado que la prueba arrojo un valor positivo se le conoce como 
Valor Predictivo Positivo (VPP) [4]. Utilizando esta simbologia 
tenemos que: 



VPP = - 



P(E)*P(R + fE) 



P(E) *P(R+fE) + P(S) *P(R + fS) 
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A la probabilidad de estar enfermo se le conoce como la 
Prevalencia [4], y un concepto muy utilizado en este campo es 
la Sensibilidad de una prueba de diagnostico. Esta representa la 
probabilidad de acertar en el diagnostico, es decir, representa la 
probabilidad de dar un resultado positivo en los enfermos [4]. De 
modo que: 

Sensibilidad = P(R + IE) 

Sustituyendo estos conceptos en la expresion del Teorema de 
Bayes tenemos: 

_ prevalencia * sensibilidad 



prexafettcki*semiMitiatt t P{S)*P{R+ t'S) 

Nos quedan todavia los terminos de P(S) y P(R+/S). 

Sabemos que E y S son dos eventos mutuamente excluyentes 
yquesecumpleque: 

P(E) + P(S) = l 

Despejando de esta expresion P(S) tenemos que: 

P(S) = 1-P(£) 6 

P(S) = 1 - prevalencia 

TASA DE LOS FALSOS POSITIVOS 

Por otro lado, P(R+/S) representa la probabilidad de dar un 
resultado positivo en las personas sanas. Esta probabilidad es 
conocida como la tasa de los falsos positivos y es un tipo de error 
que se comete en las pruebas de diagnostico [5]. 

De modo que la expresion para VPP puede expresarse como: 

Los valores de la tabla 1 pueden reescribirse en terminos de sus 
cardinalidades. Esto se muestra en la tabla 2. 

prcvaie ncia * semibtlidad 



prmtalmtda * wmibitidad + (1 - prevalencia)* tasa dc fais&s poslttv» 



Tenemos que: 



Resultado 
de la prueba 


Enfermo 

E 


Sano 
S 


R+ 


VP 


FP 


R- 


FN 


VN 



Tabla 2 Presenta los valores en un estudio de casos 
y controles. 



VP: significa verdaderos positivos, es decir, estan enfermos y la 
prueba les dio positivos. 

VN: verdaderos negativos; estan sanos y la prueba les dio 
negativo. 

FP: falsos positivos; estan sanos y la prueba les dio positivo. 

FN: falsos negativos; son pacientes enfermos que la prueba les 
da negativo. 

Otro concepto muy mencionado en las pruebas de 
diagnostico es la especificidad. Esta representa la capacidad que 
tiene la prueba para detectar la ausencia de la enfermedad en 
personas sanas [5], esto es, P(R7S), llamada tambien la tasa de 
los verdaderos negativos. Por la definicion de la probabilidad 
condicional tenemos que: 



especificidad = P(R - / S) = 



P(R-nS) 
P(S) 



Observe como la P(R+/S), que es la tasa de los falsos positivos 
y la especificidad son calculadas sobre el total marginal de 
personas sanas, de modo que sus valores son complemento. Es 
decir, 1 -especificidad nos da la tasa de los falsos positivos. 

Asi que la expresion para el valor predictivo positivo queda 
concluida como 



vpp- 



prevalencia * sensibilidad 



prevalencia * sensibilidad +■ (1 - prevalencia) *(l - especificidad) 



Observe como la probabilidad de que el paciente realmente 
tenga la enfermedad, dado que la prueba da positivo, esta en 
funcion de la prevalencia, sensibilidad y especificidad, conceptos 
conocidos por el medico. 

VALOR PREDICTIVO NEGATIVO (VPN) 

De igual forma, podemos definir la probabilidad de estar sano, 
dado que la prueba da un resultado negativo. A este concepto 
se le llama Valor Predictivo Negativo (VPN). Expresado como 
probabilidad condicional, tenemos que 



VPN = P(S/R-) 



P(SnR-) 



Observe en la tabla 3 como tambien se cumple que: 



P(R-) - P(E n R-) + P(S n R-) 



Resultado 
de la prueba 


Erfermo 

E 


Sano 
3 




R+ 


P(E^R+) 


P(S^R+) 


P(R+> 


R- 


P(E -R-) 


P(S-R-) 


P(R-) 




P(E) 


P(S) 
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Tabla 3 Presenta las probabilidades conjuntas y 
marginales 



Asi que, aplicando de nuevo el Teorema de Bayes tenemos: 



VPN: 



P(S)*P(R-fS) 



VPN 



VPN 



P(S) * P{R - / 5) + P(E) * P(R - / E) 

Expresando de nuevo la probabilidad de sanos en terminos 
de la prevalencia, y sabiendo que la P(R7S) es la especificidad, 
tenemos que: 

(1 - prevalencia) * especificidad 
(1 - prevalencia) * especificidad + prevalencia *P(R - IE) 

Anteriormente se habia hablado que el concepto de 
sensibilidad expresa la probabilidad de dar un resultado positivo 
entre los enfermos, es decir, P(R-/E). Esta probabilidad es llamada 
tambien la tasa de los verdaderos positivos, y a P(R-/E) , que 
representa la probabilidad de dar un resultado negativo entre 
los enfermos se le conoce como la tasa de los falsos negativos. 
Tenemos que los valores de ambas probabilidades son 
complementarios, de modo que P(R-/E) la podemos expresar en 
terminos de la sensibilidad, es decir, P(R-/E)= 1 -sensibilidad. Asi 
que nuestra expresion para el Valor predictivo negativo es: 

(l - prevalencia)* especificidad 



(1 - prevalencia)* especificidad + prevalencia *{1 - sensibilidad) 



Asi pues, conocer la probabilidad de que un paciente este 
enfermo dado el resultado positivo de la prueba, o no tener la 
enfermedad dado un resultado negativo de la prueba, depende 
de la prevalencia de la enfermedad asi como de la sensibilidad y 
especificidad de la prueba. 

UNA APLICACION HIPOTETICA. 

Se hizo un estudio de casos y controles, donde los casos eran 
personas que padecian de glaucoma y los controles personas 
sanas. El glaucoma es un trastorno que se presenta en el nervio 
optico y se debe al aumento de la presion en un ojo, tambien 
conocida como presion intraocular (PIO) [6]. Para medir la 
presion intraocular el medico utiliza una prueba de tonometria 
[7]. Para una persona normal, el tonometro debe registrar una 
medida entre 10 y 20 mmHg, asi que una presion intraocular 
alta se considera a partir de 21 mmHg [8]. Si el medico sabe que 
aproximadamente un 7% de la poblacion padece de glaucoma 
y que la sensibilidad de tonometro es del 87.5 y tiene una 
especificidad del 91 % <;Cual es la probabilidad de que un paciente 
realmente tenga glaucoma dado que la lectura en el tonometro 
fue mayor de 20 mmHg? 



Para este caso calculemos el VPP: 



VPP-. 



prevalencia * sensibilidad 



prevalencia * sensibilidad + (I - prevalencia) *{1 - especificidad) 



Segun la informacion proporcionada, la prevalencia de la 
enfermedad es de 0.07, por tanto, (1 -prevalencia) que es la 
probabilidad de estar sano, es de 0.93. La sensibilidad es de 
0.875 y la especificidad es de 0.91. La probabilidad de los falsos 
positivos (1 -especificidad) es igual a 0.09. Asi que la probabilidad 
de tener glaucoma dado un resultado positivo en el tonometro 
es: 



VPP-- 



(.07)* (.875) 



(.07)* (.875) + (.93)* (.09) 



r(X42 



Claro que el medico tambien debe de basar su diagnostico 
en otros factores de riesgo como son antecedentes familiares, la 
edad, raza, enfermedades como el diabetes, etc. 

Es importante mencionar que en muchos de los textos de 
epidemiologia o de bioestadistica medica expresan las formulas 
de estos conceptos de una manera mas simple. Estos se basan 
en una tabla de contingencia como las que ya se han presentado 
en este escrito, pero con una notacion diferente para los valores 
dentro de las casillas. Ver la tabla 4. 



Resultado 


Enfermo 


Sano 


Total 


R+ 


A 


B 


A+B 


R- 


C 


D 


C+D 


Totai 


A+C 


B+D 


A+B+C+D 



De modo que las expresiones para el calculo de sensibilidad, 
especificidad, valores predictivos positivos y negativos y ademas 
la prevalencia pueden expresarse como: 



sensibilidad - 



especifwidad 



VPP-- 



A + C 

D 
B + D 



A + B 



VPN 



D 
' C + D 



prevalencia = 



A + C 



A+B+C+D 



Puede verse que en las primeras cuatro expresiones esta 
implicita la probabilidad condicional. 

CONCLUSIONES 

Hasta aqui, se ha intentado presentar una aplicacion sencilla y 
sin complejas demostraciones de la probabilidad condicional y 
el Teorema de Bayes. Cabe la pena mencionar que no es la unica 
en este campo, existen otras como por ejemplo en los estudios 
de equivalencia, monitorizacion de ensayos clinicos, evaluacion 
de datos locales, etc. Asimismo, en otras areas como las ciencias 
sociales y la ingenieria, por mencionar algunas, se presenta 
tambien la aplicacion de la teoria bayesiana. 
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El CONSUMO DE ALCOHOL EN LAS ENTIDADES DEL 

NORTE DE MEXICO 
;UN PROBLEMA DE SALUD PUBLICA EN LOS ESTADOS 

FRONTERIZOS? 



SALVADOR PONCE SERRANO 



£/ consumo nocivo de alcohol en las personas tiene graves consecuencias en su salud y productividad. Se debe 
contar con cifras e indicadores para cuantificar la magnitud del problema para poder sugerir medidas correctivas 
y preventivas. Las encuestas nacionales son herramientas que contienen datos valiosos quepermiten diversos andlisis 
para caracterizar la information sobre este problema. En este documento se identifica la magnitud del problema de 
consumo de alcohol a partir de los datos de la Encuesta Nacional de Adicciones 2008, se comparan los indicadores de 
consumo excesivo en la frontera norte, en cada uno de los estados que la integran y los datos promedio del pais. La 
frecuencia de consumo excesivo encontrada en la poblacion de la frontera norte fue 60.6 hombres y 14.5 mujeres por 
cada 100 personas, en edades de 12 a 65 ahos, Aunque el consumo consuetudinario es mas tfpico de la region que en 
el resto del pais, predomina el consumo alto, sobre todo en cuatro estados, entre ellos, Sonora. 
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